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Abstract: We study how codimension-two objects like vortices back-react gravitationally 
with their environment in theories (such as 4D or higher-dimensional supergravity) where the 
bulk is described by a dilaton-Maxwell-Einstein system. We do so both in the full theory, 
for which the vortex is an explicit classical ‘fat brane’ solution, and in the effective theory 
of ‘point branes’ appropriate when the vortices are much smaller than the scales of interest 
for their back-reaction (such as the transverse Kaluza-Klein scale). We extend the standard 
Nambu-Goto description to include the physics of flux-localization wherein the ambient flux of 
the external Maxwell field becomes partially localized to the vortex, generalizing the results of 
a companion paper [4] to include dilaton-dependence for the tension and localized flux. In the 
effective theory, such flux-localization is described by the next-to-leading effective interaction, 
and the boundary conditions to which it gives rise are known to play an important role in 
how (and whether) the vortex causes supersymmetry to break in the bulk. We track how 
both tension and localized flux determine the curvature of the space-filling dimensions. Our 
calculations provide the tools required for computing how scale-breaking vortex interactions 
can stabilize the extra-dimensional size by lifting the dilaton’s flat direction. For small vortices 
we derive a simple relation between the near-vortex boundary conditions of bulk fields as a 
function of the tension and localized flux in the vortex action that provides the most efficient 
means for calculating how physical vortices mutually interact without requiring a complete 
construction of their internal structure. In passing we show why a common procedure for doing 
so using a d-function can lead to incorrect results. Our procedures generalize straightforwardly 
to general co-dimension objects. 
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1. Introduction 

A ‘vanilla’ Nielsen-Olesen vortex [1, 2] carries U{ 1) flux but this flux is normally expelled 
from the region outside of the vortex. In this paper we study the gravitational response of 
vortices (or ‘fat’ branes) that carry localized amounts of an external magnetic flux that is not 
expelled from its surrounding environment (so-called ‘Dark Vortices’ or ‘Dark Strings’ [3]). 
Our description of these systems closely parallels our companion paper [4], extending it to 
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the case where effective couplings are functions of the bulk dilaton that tends to appear in 
supersymmetric theories. 

Vortices which partially localize bulk flux can arise within supersymmetric theories in 
various dimensions, and whether or not their presence breaks supersymmetry depends on the 
relative size of their tension and the amount localized bulk flux they carry [5]. Because of 
this their tension and localized flux compete with one another in the amount of curvature 
produced by their back-reaction on their surrounding geometry. Our main goal is to explore 
this competition in detail and to identify precisely how it depends on the various parameters 
that describe the vortex physics. 

We have several purposes in mind when doing so. First and foremost we wish to under¬ 
stand how brane back-reaction influences the transverse geometry through which the vortices 
move, and in particular how they stabilize the value of the dilaton and so set the size of the 
transverse dimensions and the curvature of the geometry induced on their world-sheet. Much 
is known about the systems in the limit when the sources are pointlike [6] , and in particular it 
is known that the world-sheet geometries are exactly flat (at the classical level) if the dilaton 
should have a vanishing derivative at all source positions [7] (a result which we also re-derive 
here). 1 What this leaves open is whether there exists any kind of source for which a vanishing 
near-source dilaton derivative is possible and, even if so, whether the resulting curvature is 
then nonzero but dominated by finite-size vortex effects that could be suppressed for small 
vortices but not vanishing. 

We find three main results. 

• Vortex-dilaton coupling: In general the radial derivative of any bulk field very near a 
point source is completely controlled by the derivative of the source action with respect 
to the field of interest [6, 8, 9] (much like the quantity lim r _»o r 2 dV/dr is dictated by 
a point-source’s charge, Q oc 5S/SV in electrostatics). So naively a vanishing dilaton 
derivative at a the position of a source brane is arranged by not coupling the dilaton to 
the brane at all. While we confirm the truth of this assertion, we also find it is harder 
to completely avoid such a brane coupling to the dilaton than one might think. More 
specifically, we find that although it is possible to arrange a dilaton-free tension, it is 
much more difficult to arrange dilaton-free localized flux. It is more difficult because in 
a supersymmetric bulk the value of the dilaton sets the local size of the gauge coupling 
for the flux. 

• Modulus stabilization: By computing how vortex-bulk energetics depend on the value of 
the dilaton we verify earlier claims [6, 8, 10] that (with two transverse dimensions) brane 
couplings generically stabilize the size of the transverse dimensions in supersymmetric 
models, in a manner similar to Goldberger-Wise stabilization [11] in 5D. They do so 
because they break the classical scale invariance of the bulk supergravity that prevents 

1 Indeed this underlies the study of these system as potential approaches [12, 13] to the cosmological constant 

problem [13, 14]. 
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the bulk from stabilizing on its own (through eg flux stabilization). The tools we provide 
allow an explicit calculation of the energetics as a function of the dilaton held (and so 
in principle allow a calculation of the stabilizing dilaton potential). 

• Low-energy on-brane curvature: We find that the same dynamics usually also curves 
the dimensions along the vortex world-sheets, and generically does so by an amount 
commensurate with their tension, R ~ G N T , where T is the vortex tension (defined 
more precisely below) and G N is Newton’s constant for observers living on the brane. 
For specific parameter regimes the on-vortex curvature can be less than this however, 
being suppressed by the deviation of the vortex from scale invariance (when this is 
small) and the ratio of the vortex size to the 

• Matching and effective descriptions: We describe our analysis throughout in two com¬ 
plimentary ways. On one hand we do so using the full (UV) description within which 
the vortices are explicit classical solutions. We then do so again using the lower-energy 
(IR) extra-dimensional effective theory within which the vortices are regarded as point 
sources because their sizes are not resolved. By comparing these calculations we pro¬ 
vide explicit matching formulae that relate dilaton-dependent effective parameters (like 
tension and localized flux) to underlying properties of the UV completion. 

• Efficient description of point-source back-reaction: We provide explicit formula that re¬ 
late the near-source boundary conditions of bulk fields in terms of the source tension and 
localized flux. Because these boundary conditions determine the integration constants 
of the external bulk solutions they efficiently solve the back-reaction problem in a way 
that does not depend on providing a detailed construction of the internal microstructure 
of the sources. As such they provide the most efficient way to describe how such sources 
gravitate, and the framework within which to renormalize the divergences associated 
with the singularity of bulk fields at the source positions [10, 15]. We show in passing 
why a commonly used (5-function way of trying to infer brane-bulk interactions can give 
incorrect results. 

Our conclusions also include several more technical observations about the gravitational 
physics of small localized brane sources. In particular, we identify how stress-energy conser¬ 
vation strongly constrains the components of the source source stress energy tensor, allow¬ 
ing in particular the extra-dimensional off-brane stress-energy components to be dictated in 
terms of the effective action governing the on-brane degrees of freedom. We do so by using 
the vortex system to explicitly construct the effective theory of point sources describing the 
extra-dimensional response to the vortices on scales too long to resolve the vortex size. In 
particular we show how the radial Einstein ‘constraint’ determines the two nontrivial off-brane 
components of stress energy (for rotationally invariant - monopole - sources) on scales much 
longer than the size of the vortices themselves. 
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A road map 

We organize our discussion as follows. The following section, §2, describes the bulk system 
in the absence of any localized sources. The purposes of doing so is to show how properties 
of the bulk physics (such as extra-dimensional size and on-brane curvature) are constrained 
by the field equations, which controls the extent to which they depend on the properties of 
any source branes. 

Then come sources. §3 first describes the source physics in terms of a UV theory within 
which the branes can be found as explicit classical vortex solutions. This is followed, in §4 by 
a discussion of the higher-dimensional effective theory that applies on length scales too large 
to resolve the vortices, but small enough to describe the extra-dimensional geometry. In this 
regime the vortices are described by effective point sources, and we find their properties as 
functions of the choices made in the full UV theory. §5 summarizes how the main physical 
quantities scale as functions of the couplings assumed between the vortex and the dilaton. 
Our conclusions are summarized in a final discussion section, §6. 


2. The Bulk 


We start by outlining the action and field equations of the bulk, which we take to be an 
Einstein-Maxwell-scalar system. Our goal is to understand how bulk properties (such as 
curvatures) are related to the asymptotic behaviour of the fields near any localized sources. 
We return in later sections describing the local sources in terms of vortices, with the goal of 
understanding what features of the source control the near-source asymptotics. We imagine 
the bulk to span D = d + 2 spacetime dimensions with the d-dimensional sources localized 
in two transverse dimensions. The most interesting cases of practical interest are the cosmic 
string [with ( D,d ) = (4,2)] and the brane-world picture [with ( D,d ) = (6,4)]. Since this 
section on the bulk duplicates standard results [7, 12, 16, 17], aficionados in a hurry should 
skip it and go directly to the next section. 


2.1 Action and field equations 

The bulk action of interest is given by 


S B = - d d+2 x 


7^2 g MN (tZ-mn + + V B {(f>) + -j-e ^A mn A a 


J d d+2 x \J—g ^ L eh + + L/^j —: — j d d + 2 x y/—g L f 


( 2 . 1 ) 


where 2 A MN = d M A N — d N A M is a D-dimensional gauge Held strength and 1Z MN denotes the 
D-dimensional Ricci tensor. The second line defines the Einstein-Hilbert, scalar and gauge 

2 We use Weinberg’s curvature conventions [18], which differ from those of MTW [19] only by an overall 
sign in the definition of the Riemann tensor. 
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contributions — i.e. L EH , L^ and L A — in terms of the items in the line above, with L B 
denoting their sum. When needed explicitly we take the scalar potential to be 

V b (4>) = Vo , (2.2) 

and in the special case Vo = 2g\/ k 4 this corresponds to a subset of the action for Nishino- 
Sezgin supergravity [20] when d = 4 (so D = 6). In this case g R is the gauge coupling constant 
for a specific gauged U(1) R symmetry that does not commute with 6D supersymmetry. In 
what follows we denote by g A the gauge coupling for the gauge held A M , although g A = g R 
in the most interesting 3 situation where this gauge held is the U(1) R gauge held. 

Scaling properties 

For later purposes notice that L B scales homogeneously, L B —> s~ 1 L B under the rigid rescal¬ 
ings g MN —> s g MN and —> s _1 e^, which as a consequence is a symmetry of the classical 
equations of motion. The corresponding transformation for the action is S B —> s d ^ 2 S B , as 
is most easily seen by transforming to a scale-invariant metric g MN = e^g MN , in which case 
all terms of the bulk lagrangian are proportional to e~ d< bl 2 with cf> otherwise only appearing 
through its derivative, d M <p. Besides ensuring classical scale invariance this also shows that 
it is the quantity e d< that plays the role of h in counting loops within the bulk part of the 
theory. 

The bulk system enjoys a second useful scaling property. If we rescale the gauge field 
so A m : = g A A M then arbitrary constant shifts 0 —> (j) + leave the action unchanged 
provided we also rescale the couplings by g A — > g Ait := g A e and Vo — > 14 := Vo e^* (or, if 
Vo = 2g B /K, 4 , equivalently g R —> g R * := g R e This is convenient inasmuch as <fi = 0 can 
always be chosen to be the present-day vacuum provided the values of constants like g A and 
Vq are chosen appropriately. 


Bulk field equations 

The field equations for the Maxwell field arising from the bulk action are 

9m (V-9 e"^ MJV ) = 0 , (2.3) 

which is supplemented (as usual) by the Bianchi identity d A^) = 0, for the 2-form A MN . The 
bulk dilaton equation is similarly 


= -J= g d M g MN d N (k) = k 2 (v B - l a ) , 

while the Einstein equations can be written in their trace-reversed form 


(2.4) 


TZ-mn — K Xt 


(2.5) 


3 Enhanced interest comes from the unbroken 4D supersymmetry that these configurations can enjoy. 
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where X MN := T MN — (1/d) g MN T p P and the bulk stress-energy tensor is 


rp{ B ) 
- L MN 


—o ^A mp A n p — gMN 

i/.Z 


'(W 

2k 2 


+ V B + Ta 


( 2 . 6 ) 


Notice that it is a special feature of the split into D = d + 2 dimensions that any maximally 
symmetric contribution to the d-dimensional part of the stress-energy, oc g^, drops out 
of and so naively does not contribute directly to the 4D Ricci curvature, TZ^. 


2.2 Symmetry ansatze 

Our interest throughout is in solutions that are maximally symmetric in the d ‘spacetime’ (or 
‘on-brane’) dimensions (spanned by x M ) and axially symmetric in the 2 ‘transverse’ dimensions 
(spanned by y m ). We therefore specialize to fields that depend only on the proper distance, 
p, from the points of axial symmetry, and assume the only nonzero components of the gauge 
field strength, A mn , lie in the transverse two directions. We choose the metric to be of the 
general warped-product form 

ds 2 = g MN dx M dx N = g mn dy m dy n + g, w dx^dx u , (2.7) 

with 

9mn = 9mn{y) and 9n V (x,y) = W 2 {y)g llv {x ), (2.8) 

where g ia ,(x) is the maximally symmetric metric on d-dimensional de Sitter, Minkowski 
or anti-de Sitter space. The corresponding Ricci tensor is 1Z MN dx M dx N = 1Z^ dx^dx v + 
'JZmn dy m dy n , and is related to the Ricci curvatures, and R m n , of the metrics g^ and 
9mn by 

TZ^u = R„ v + g mn 

and 

TZmn = Rmn + V m V n W , (2-10) 

where V is the 2D covariant derivative built from g mn . For the axially symmetric 2D metrics 
of interest we make the coordinate choice 


(d - 1 )d m Wd n W + WV m V n W 


9^v > 


(2.9) 


gmn dy m dy n = A 2 {r) dr 2 + B 2 {r) d 9 2 = dp 2 + B 2 {p) d0 2 , (2.11) 

where proper distance, p, is defined by dp = A(r)dr. Some useful properties of these geome¬ 
tries are given in [4], 

With these choices the held equation simplify to coupled nonlinear ordinary differential 
equations. Denoting differentiation with respect to proper distance, p, by primes, the gauge 
held equations become 

(e-*W d A p o 



B 


( 2 . 12 ) 








while the field equation for the dilaton becomes 

(BW d <!>')'= k 2 (v b -L a ) =:k 2 X b , (2.13) 

where the last equality defines X B := V B — L A . With the assumed symmetries the nontrivial 
components of the matter stress-energy are 


Tf_ w = —gftu Qb , T p p = Z — X and T e e = -{Z + X), (2.14) 

where the bulk contribution to X is X B as defined above and 

q b := ( ^- + Vb + La and Z B . (2.15) 

The three independent components of the trace-reversed bulk Einstein equations then 
consist of those in the directions of the d-dimensional on-brane geometry, 

'R'fiv = ^ k 2 X b g pv , (2.16) 

of which maximal symmetry implies the only nontrivial combination is the trace 


R 


R-(d) ■— g^R'nv — w2 + 


d 


W 2 BW d 


BW'W 


d -1 


= -2 k 2 X R 


(2.17) 


and we use the explicit expression for 7 Z^) i n terms of R and W. 

The components dictating the 2-dinrensional transverse geometry similarly are 7 Z mn = 
—K 2 X mn , which for the bulk has the following two nontrivial components: 


:w" b’W 

Km ■■= g mn n mn = R+d(— + w - 


= -K 2 X m m = -2 k 2 


Qb - [1 — -) X B 


, (2.18) 


and the other can be the difference between its two diagonal elements 


Q% - G\ = n p P - n\ = -« 2 ( t p p - t%) , 


which for the assumed geometry becomes 


B fW r \ 
W \lf J 



Ml 2 

d 


< 0 . 


(2.19) 


( 2 . 20 ) 


This shows that W'/B is a monotonically decreasing function of p. 

A useful linear combination of the above Einstein equations corresponds to the (89) 
component of the trace-reversed equation, which for the bulk reads 


(B'w d y 

BW d 
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2.3 Bulk solutions 


We next describe some of the properties of the solutions to these field equations. In order to 
accommodate our later inclusion of the equations governing any localized sources, we examine 
solutions of the bulk equations only within a domain, the ‘bulk’: B ex t, which consists of the 
full 2D geometry transverse to the sources from which small volumes (‘Gaussian pillboxes’, 
X v ) are excised. These pillboxes completely enclose any sources that might be present. In 
practice we define B ext such that the radial proper distance coordinate satisfies p- < p < p+, 
with the sources lying just outside of this range. When not specifying which source is of 
interest, we generically use p v = {p +1 p-} to denote the boundary between the source and 
the bulk. 

Integral relations 

Before writing some exact and approximate solutions we first record several exact integral 
expressions that can be derived by directly integrating the field equations over the volume 
B ext , being careful to keep track of its boundaries at p = p- and p = p + . 

Integrating the bulk Maxwell equation, (2.12), with respect to p in radial gauge gives 

A p e -A e - , (2.22) 

for integration constant Q, and this in turn integrates locally to give (up to a gauge transfor¬ 
mation) 

Mp + ) - Mp-) = £ • < 2 - 23 ) 

where we define the notation 

P+ 

(" ' L : = (...)= 2, / dp (...) . (2.24) 

P — 

These expressions also contain the seeds of flux quantization when applied to spherical 
transverse dimensions. To see this we take p± to he infinitesimally close to the north and south 
poles and excise these two points, leaving the topology of a sphere with two points removed 
(or an annulus). Integrating (2.23) around the axial direction and using the quantization of 
any gauge transformation, g~ l deg, then implies 

77 = / Ae df) - j A„ at) + Q ( = 4 >,_ + + Q ( ^ 53 )^ , (2.25) 

P- P+ 

where N is an integer and g A is the gauge coupling for the held A M and 3> A ± denotes the 
total A flux through the caps over the relevant poles. In the absence of sources at the poles 
we can contract the circles at p = p± to a point and learn j> ± Aq d0 — > 0, in which case (2.25) 
becomes the usual condition on Q required by flux quantization. But when sources are present 








<f ± Aq dd need not vanish if there is flux localized within the source. When this is so their 
presence on the right-hand side of (2.25) modifies the implications for Q of flux-quantization 
[ 6 ], 

Similarly integrating the field equation (2.13) for the dilaton gives 


BW d ft 

— 

BW d (j)' 

- 

P+ 

■ 


J p- 


n 2 k 2 

(y^ B ^' 4 )ext ‘ 


(2.26) 


Two of the Einstein equations also provide direct first integrals. Integrating (2.21) leads to 

.2 / / \ - 2 


B'W d 


j p+ 


B'W d 


K 


J P- 


2vr 


= -—(q b -Z b -[1--X j 


d 


= --(l* + w 

ext \ / ext 


(2.27) 


while the integral of (2.17) implies 


1 =-i[fl<M , T.„ + 2 '‘ 2 <-r»L,] ■ (2.28) 

J p + L J p_ 


Of particular interest for present purposes is the simple relationship between the on-brane 
curvature, R, and the near-source boundary values of bulk fields obtained by combining (2.26) 
with (2.28): 


BW d 




BW d 





ext 


(2.29) 


This states that the quantity BW d (^ + f In W) 1 is monotonic in p everywhere outside of 
the sources, growing or shrinking according to the sign of R or remaining constant when 
is flat. Should BW d (</> + | In IT) take the same value for two different values of p in 
the bulk, then we can conclude that R = 0. In particular, if there exist sources for which 
2 (f>' + dW'/W = 0 at both source positions, then any interpolating geometry between the 
sources must satisfy R = 0; that is, the vanishing of (e^W d ^ 2 ) in the near-source limit for 
both sources is a sufficient condition for R = 0 (as first argued some time ago [7]). A special 
case of the R = 0 solutions are those for which BW d (</> + | In IT) vanishes everywhere, for 
which W d / 2 is constant. 


Asymptotic forms 

Given the intimate relationship between R and near-source derivatives implied by eq. (2.29), 
we next examine what the field equations imply about the form of the bulk solutions very 
close to, but outside of, a source much smaller than the transverse dimensional size (eg outside 
a source for which p v —> 0). 4 

4 Our treatment here follows closely that of [4] and [21], which are themselves based on the classic BKL 
treatment of near-singularity time-dependence [22]. 
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Usually the presence of such a small source induces an apparent singularity into the 
external geometry, at the position p = pj’ defined as the place where B(p) vanishes when ex¬ 
trapolated using only the bulk held equations. That is, it implies there would be a singularity 
in the external geometry if this geometry were extrapolated right down to zero size (rather 
than being smoothed out by the physics of the source interior, as illustrated in Figure 1). If 
the centre of the source is chosen to be p = 0 then in order of magnitude p* is of order the 
source’s size: p* ~ p v . 



Figure 1: A cartoon illustration of the definition of p*. The (blue) metric function B increases 
linearly away from the origin with unit slope B(p) ~ p. Outside of the source p > p v , the solution is 
a power law in p with B(p) ~ p b . The straight (red) line extrapolates this exterior behaviour to the 
point, p = p*, where the external B would have vanished if the vortex had not intervened first. 


When constructing a near-source solution it is most informative to do so as a series 
solution expanding in powers of the distance, p := p — p*, from the singular source point. 
That is, 


\ W / Ul + l 

V +w iv + - 



(2.30) 


where l is a measure of the proper size of the transverse geometry, which we assume to be 
much larger than the source’s size, so £ p. The powers w, b and z describe the nature of 
the singularity at p = p*, and all but three combinations of these parameters and the W t , 
Bj and 4>i coefficients turn out to be related to one another by the bulk field equations. The 
three ‘free’ parameters are instead determined by boundary conditions that the bulk solutions 
satisfy in the near-source regime. 

5 Although p* does not carry a source index v, it is clear that there is one such singular point for each source 
labelled by v. 
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In particular, for small sources these field equations allow all of the Wj, Bi and <+ to 
be computed in terms of (say) 4>o, Wo and Bo , and further imply the following two ‘Kasner’ 
relations [21, 23] amongst the powers 6, w and 2 : 

diu + b = dw 2 + b 2 + z 2 = 1. (2-31) 


The second of these in turn implies w, b and 2 must reside within the intervals 

|u;| < —j= and |&|,|z|<l. (2.32) 

V d 

It turns out that the rest of the field equations do not give additional constraints on the three 
parameters b, w and z, and so one combination of these is one of the quantities determined 
by the physical properties of the source. 

Expansion about a regular point — ie not the location of a singularity — should corre¬ 
spond to the specific solution z = w = 0 and b = 1 to eqs. (2.31). In the presence of weakly 
gravitating sources we expect to find small deviations from these values, whose size can be 
inferred by solving (2.31) perturbatively. To this end we write 6=1 + b\5 + 62 < 5 2 + • ■ • , 
w = w\5 + • • • and z = z\5 + • • • and expand. Working to order 5 2 we find the one-parameter 
family of solutions 

2 2 

z = z 1 S + 0(S 2 ), 6 = 1 - Z — + 0{5 3 ) and dw = y + 0(5 3 ). (2.33) 

Quite generally only 2 can deviate from its background value at linear order, and the leading 
quadratic contributions to w and 6—1 are not independent of this linear deviation in 2 . In 
later sections we find 2 : is determined by the coupling of cf) to the source lagrangian, and so 
we generically expect 5 to be of order k 2 v 2 , where v 2 is a measure of the energy density (or 
tension) carried by the source (more about which below). Fig. 2 provides numerical evidence 
for the validity of the Kasner equations and their perturbative solution, (2.33), nearby but 
outside of an explicit vortex solution. 

In terms of these Kasner exponents, the combination appearing on the left-hand side of 
eq. (2.29) in the near-source limit is 


lirn BW d 

p^p* 



(2z + dw) 


B 0 W 0 d 

£ 



dw-j-b— 1 


(2z + dw) 


BpW d 

£ 


(2.34) 


where the second equality uses the linear Kasner condition (2.31). When eq. (2.34) is order 
unity then eq. (2.29) implies that the curvature R is of order R ~ (k^/k 2 ) corresponding to a d- 
dimensional energy density of order 1 /k, 2 . Here Kd is the dimensionally reduced gravitational 
coupling of the low-energy d-dimensional theory, and we use a result, derived below, that 
n~i 2 - « -2 (W^ 2 )ext- 

Of particular interest are situations where the near-source solutions satisfy 2 ; = 0 since the 
Kasner conditions imply these also satisfy w = 0 and so 22 : + dw = 0. Consequently for any 
such source the leading contribution, eq. (2.34), to R vanishes. In this case it is a subleading 
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In(pAv) 


In (p/r v ) 


Figure 2: This plot illustrates the near-source Kasner region for bulk fields that are sourced by 
a generic vortex with tension v 2 = 0.25/re 2 . In the top left plot, the vortex profiles are shown on 
a logarithmic axis, with the (blue) scalar profile F increasing towards its asymptotic vacuum value 
F —>• 1 and the (grey) gauge profile P decreasing towards its asymptotic value P —> 0. (See below 
for more detailed definitions.) These profiles approach their vacuum values exponentially, and they 
can be neglected in the region r = ln(p/r v ) > 2, which is where the Kasner region p > p v begins. In 
the next plot, ln(B') is seen to be linear in r in the darkened (blue) near-vortex region, before the 
bulk sources dominate the behaviour of B' at larger r. This blue region is the Kasner region, where 
the bulk fields obey power laws, and this behaviour is also apparent in the plots of ln(e^) = <f> and 
ln(IT d ) = din IT. The Kasner powers can be extracted in this region, and for this particular source 
we find numerically dw = 0.0034, 1 — b = 0.0034 and z = 0.082 which is in agreement with the 
perturbative solution of (2.33). Finally, we note that these numerical results are also consistent with 
the estimates z ~ S ~ re 2 ?; 2 as is argued below. 


term that first contributes on the left-hand side of (2.29), implying that R is suppressed 
relative to 1/re 2 by a power of the ratio of source and extra-dimensional sizes: p±/&. This 
asymptotic reasoning is borne out by the explicit numerical and analytic solutions described 
in the next sections. 

Explicit bulk solutions 

It is useful to see how the above general arguments go through for explicit solutions to the 
bulk field equations, which are known in great detail [16, 17] when <// —y 0 in the near-source 
limit in the special case where d = 4 and D = 6 and Vq = 2g\/n A . (See also [24] for a 
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discussion of linearized solutions; [21, 25] for solutions with 6 (p 1 ^ 0; [28] with other fields 
nonzero; [29] for time-dependent configurations; [30] with more than two brane sources; and 
with black hole solutions [31].) 

The solutions are most simply written using the symmetry ansatz 

ds 2 = W 2 (£) ds 2 4 + r 2 (£) (d£ 2 + a 2 (£) sin 2 £ df? 2 ) , (2.35) 


where ds 2 denotes the maximally symmetric on-vortex geometry, ds 2 = g^ v dx^dx". With 
this ansatz, as seen above, the field equations ensure R = 0 provided we assume <$ —> 0 in 
the near-source limit, which we now do (and so also take g^ v = g^u)- The dilaton and metric 
function then turn out to be 


e m 




o(sn w\i) 


and r(£)=r B W(£)e ^ 2 , 


(2.36) 


with 


IT 4 (£) = e v sin 2 ^ + e v cos 2 ^ = cosh v — sinh v cos £ 

where v, T and p are integration constants. Notice that r 2 e^ = r| 
length-scale r B set by 


K 



, (2.37) 

for all £, with the 

(2.38) 


The two sources for this geometry are located at the two singular points, £_ = 0 and 
£ + = 7 r, where the transverse space pinches off. In the near-brane limit, £ —> 0, we have 
W —> W- + C7(£ 2 ), and so the proper radial distance, given by 


P(£) = J d£ r(£) = r B e~^ 2 f d£ W(£), 


(2.39) 


in this limit becomes p = r B e - ^/ 2 [Wo£ + 0(£ 3 )], and so a — > a- = T/Wf + 0(p 2 ), B — > 
a-p + 0{p 3 ) and e^ /W 2 + 0(p 2 ) in the near-brane limit (corresponding to the Kasner 

exponents z = w = 0 and 6 = 1 ). Similar properties also hold for £ — > tt. 

Two of the integration constants — v and T — can be traded for the conical defect 
angles, <5± = 27r(l — a±), in the two near-brane limits, with 


T = yTaqTaTr 



(2.40) 


In terms of these a(£) and W(£) are given simply by 

-j-r = — cos 2 - + — sin 2 - and IT 4 (£) = IT 4 cos 2 - + Wl sin 2 - . (2-41) 

a(£) a- 2 a + 2 2 2 

6 The de Sitter solutions to these equations are interesting in their own right as a counter-examples [26] to 
no-go theorems for the existence of de Sitter solutions in supergravity [27]. 
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In particular, in the special case W + = W- the function W(£) — and so also <?)(£), r(£) and 
a(£) — becomes constant, and the geometry reduces to the simple rugby-ball solution [12]. 
This has proper distance p = £r B eT^! 2 and 


B(p) = ar B e sin 




(2.42) 


and so £ = r B e - ^/ 2 = ^n/g R (ip) — where g R {g>) = g R e' f ^ 2 — represents the proper ‘radius’ 
of the transverse dimensions, whose physical volume is fl = A-kI 2 . 

The gauge field for the general solution with different brane properties is given by 


_ QT sin£ 
2g A W*(Q 


(2.43) 


where g A is the corresponding gauge coupling constant. Comparing with (2.22) shows Q is 
related to Q by Q = Q/(2g A r 2 B ) and so is not an independent constant. The total amount of 
flux present in the bulk (region B ext ) then is 


/ 


-4(2) 



d£dd A& 


2ttQT 

9a 


(2.44) 


so flux quantization, (2.25), relates any source flux to the defect angles by 

, r _ 9a^a- 9a^a+ gy, 

2tt 2tt ^ 


(2.45) 


This shows why brane-localized flux is generically necessary to satisfy flux quantization for 
two branes with generic tensions if Q is otherwise fixed. When there is no localized flux, we 


can estimate 


Q ~ 


1 

9aT I 


r^j 



rs j 


9r 

K 2 


(2.46) 


3. Sources - the UV picture 

The previous section describes several exact consequences of the bulk field equations that 
relate bulk properties to the asymptotic near-source behaviour of various combinations of 
bulk fields. In particular it shows how the on-source curvature, R, is determined in this 
way purely by the near-source combination BW d (0 + | lnIT) / , and so vanishes in particular 
when (j)' and W' approach zero in this limit. 

This section now turns to the question of how these derivatives are related to source 
properties, extending the results of [4] to include dilaton couplings and extending those of 
[32] to include nonzero brane-localized flux (more about which below). In this section this 
is done by making an explicit construction of the sources within a UV completion, as a 
generalization of Nielsen-Olesen vortices [1, 2], We do so by adding new scalar and gauge 
fields that admit such vortex solutions, with a view to understanding in more detail how near¬ 
source behaviour is controlled by the source properties. Because our focus here is mostly on 
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classical issues we do not explicitly embed the new sector into a supersymmetric framework, 
but we return to this issue when considering quantum corrections in subsequent analysis [33]. 

A key assumption in our discussion is that the typical transverse vortex size, f v , is much 
smaller than the size, I, of the transverse external space: r v <C t. Subsequent sections then 
re-interpret the results found here in terms of the ID-dimensional IR effective theory applying 
over length scales f v <C p < l. We follow closely the discussion of [4], highlighting the 
differences that arise as we proceed (the main one of which is the presence of the dilaton zero 
mode). 


3.1 Action and field equations 

We start with the action and field equations for the UV completed system describing the 
sources. With Nielsen-Olesen solutions in mind, we take this ‘vortex’ — or ‘brane’ — sector 
to consist of an additional complex scalar coupled to a second 1/(1) gauge field. Again we 
work in D = d + 2 spacetime dimensions, with the cases ( D,d ) = (4,2) and ( D,d ) = (6,4) 
being of most later interest. 

The full action now is S = S B + S v with S B as given in (2.1) and the vortex part of the 
action given explicitly by 

S v = - j & d+2 x y/—g 

=: ~ J < ^ Jr ~ x V~9 (j-'z + Tmix + Lq, + VbJ =: — J & dJr2 x \J—g L v , (3-1) 

where D M ^> := d M 'k — and the second line defines the various L*. The terms L z , 

Lq, and Vb describe scalar electrodynamics and are chosen to allow vortex solutions for which 
the Z MN gauge flux is localized. The L m ; x term is chosen — following [4] (see also [35, 3]) 
— to kinetically mix Z with the bulk gauge field [36] and thereby generate a vortex-localized 
component to the exterior A MN gauge flux. 

We follow [4] and write \/2 T = and adopt a unitary gauge for which the phase, fl, 
is set to zero, though this gauge will prove to be singular at the origin of the vortex solutions 
we examine later. In this gauge the term Lq, in S v can be written 


Lq, = 1 (d M il)d M il) + eVZ M Z M ) 


(3.2) 


and the vortex potential becomes 

Vbifatl’) = - -y 2 ) • 


(3.3) 


Our interest in what follows is largely in how the dependence on cj) in these interactions back- 
reacts onto the properties of the bulk, and affects the interactions of the various low-energy 
effective descriptions. We notice in passing that S v shares the classical scale invariance of S B 
only in the special case: p = r = — 1 and q = 1. 
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It is also useful to group the terms in the vortex and bulk lagrangians together according 
to how many metric factors and derivatives appear, with 


1 


Akin := x g MN 9 d m 4> d N ct> + d M ip d N tp , L gge := L A + L z + L 


1 


L vot :=V B {(t>) + V b (ct>^) and L gm :=-e 2 ^ 2 g MN Z M Z N , (3.4) 


so + Lq, + Vb = Akin + ipot + hgm- Notice that for the configurations of later interest we 
have L z > 0, L ^ > 0 and 14 > 0 while L m j x can have either sign. 


Gauge field equations 

With these choices the field equations for the two Maxwell fields are 

1 




d M 


— ( e -<t>A MN + £e r ^Z MN 


= 0, 


and 


1 




<9n. 


-g(e p *Z MN +ee r *A A 


= e 2 iJj 2 Z N , 


(3.5) 


(3.6) 


and (as usual) these are supplemented by the Bianchi identities d A = d Z = 0, for the 2-forms 
A mn and Z MN . For later purposes it is useful to write (3.5) as d M (^/^g e~ < f > A MN ^j = 0 with 
A mn defined by 

A M n '•= A mn + e e^ r+1 ^Z A 

in which case (3.6) becomes 
1 






with 


g A(0) Z MN + e(r + 1) e r *d M (j) A MN = e 2 ip 2 Z> 


A(0) ;= e P<t> _ e 2 e (2r+l)0 _ 


(3.7) 

(3.8) 


(3.9) 


We see below that the energy density of the system is given by a particular sum of the 
L^s, and when assessing the sign of the energy it is useful to notice that the off-diagonal 
contribution to L gge vanishes when this is expressed in terms of A MN rather than A MN , since 


Lgge — La + L z + L m ; x — La T L z , 


(3.10) 


where 

L a -.= -^A mn A mn and L z := ^ A(<£) Z MN Z MN . (3.11) 

This shows that the kinetic energy of the Z M gauge field is renormalized ‘ by the mixing of the 
two gauge fields, with the result only bounded below for all 7 8 real <p if p = 2r +1 and e 2 < 1. It 

7 This provides a UV perspective to what becomes a divergent renormalization [8, 9, 15] in the limit of 
zero-size sources. 

8 For some applications requiring boundedness for all cp may be too strong a criterion, since the semiclassical 
approximation relies on the assumption <C 1. Since our inference about the boundedness of the energy is 
itself performed semiclassically, it is also suspect if the unboundedness occurs only when ^ 1. If this weaker 
criterion is our guide then we only require 2r + 1 > p rather than strict equality. 
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also suggests a better split between the bulk and the vortices is to write L B + L V = L B + L\ oc , 
where L B = L EH + L ( j > + L A and L\ oc = L z + L ^ + 14 . Split this way all of the localized energy 
falls within L\ oc , because of the absence of mixing terms [4]. 

Although the quantities L A and L z are useful when describing the energy density, unlike 
in [4] their use directly in the lagrangian can lead to errors. It is important in this regard to 
keep in mind that there are two important ways in which the transition from A MN to A MN 
differs in the present case from the discussion in [4], First, if A MN is used instead of A MN in 
the field equations then one must remember that the Bianchi identity for A and Z implies A 
also satisfies 

di = e(r + l)e (r+1 ) ^d(MZ, (3.12) 

which need not vanish because of the presence in (3.7) of the field </>. The second difference 
is related to the first: one must never use eq. (3.7) to trade A MN for A MN in the action and 
then compute the field equations. This is because (3.7) is not a change of variables in the 
path integral since it is not a redefinition of the gauge potentials. In practice this kind of 
substitution is most dangerous in the <fi field equation, as may be seen from the functional 
chain rule, 


5S 

Scj)(x) 


A fixed 


5S 

A fixed 

5S 


+ I d D y 


5S 


3A MN (y) 


cj) fixed 


SA MN (y) 
6cj)(x ) 


A fixed 


■J4) +E(r + l)e^Z u J 

S n x >J Afixed V 


ss 


5 A mn (x) 


(3.13) 


infixed 


The second term on the right-hand side of this relation need not vanish when the field equa¬ 
tions are satisfied. 

For configurations with the symmetries of interest the gauge field equations reduce to 


( e-*W d A pd ^ 


= 0 , 


(3.14) 


and 


BW d 


m 


W d Z' e 

B 


+ e{r + 1 ) e 


rep 




pO 


B 2 


e 2 P 2 Z e 

B 2 ’ 


(3.15) 


where (as before) primes denote differentiation with respect to proper distance, p, and A(0) 
and A mn are as defined in eqs. (3.9) and (3.7), respectively. 


Flux quantization 


The solution to eq. (3.14) is 



QBe* 
W d ’ 


(3.16) 


which shows that A( 2 ) describes the part of the gauge fields that does not see the vortex 
sources. Ultimately the integration constant Q is fixed by the flux quantization conditions 
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for Ar 2 ) and Z( 2 ), which state 

■■= J A m = J d 2 y A p e = ^ , (3.17) 

and 

*.==/* W = /<f»»*, = -£^ = (3-18) 

where N and n tot = )T) V n v = n+ + n- are integers while e and g A are the relevant gauge 
couplings. Strictly speaking flux quantization only ensures the sum over all vortices, ntot = 
n+ + n_, is an integer. However we imagine here that the two vortices are situated at opposite 
ends of the (relatively) very large extra dimensions and so are very well-separated. Conse¬ 
quently in practice each of n + and n_ are separately integers, up to exponential accuracy. 9 

On one hand, for 4> ~ (f> v approximately constant across the narrow width of each vortex, 
this implies 



while on the other hand the left-hand side is related to Q by 

$4 = <2 J d 2 y (^ 4 ) =Q&- 4 , (3.20) 

where we define the useful notation 

:= j d 2 y^W k e* = J d 2 y ^f 2 W k , (3.21) 

to represent the 2D integrals that arise here and in later calculations. Here Q k is the integral 
of W k over the transverse dimensions using the scale-invariant metric, g mn := g mn , and 
the particular case k = 0 represents the extra-dimensional volume, H := Ho, as measured by 
this metric. 

We see Q is given by 


Q = 





£ 

e 


E 


n v e 


(r+l)<l> v 


(3.22) 


In the special case where (j) = 
becomes 

Q 


<p takes the same value at all of the vortex positions this 

2 n 
H_4 


— — ( £ ntot gC+iV 
9 a \ e ) 


(3.23) 


9 This could also be alternatively arranged by having two copies of the vortex sector, with each vortex 
carrying flux from a different 17(1) (which would therefore be separately quantized) but with both mixing with 
the bulk U( 1). 
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Scalar field equations 

The vortex scalar field equation in unitary gauge becomes 

d M (y~9 g MN d N ijj S ) = e 2 i’Z M Z M + A (y 2 - v 2 
while the dilaton equation is 

Ocj) = —L= d M (y^g g MN d N (^j = k 2 (v B - L A + qV b + pL z + rL mix ) 


(3.24) 


(3.25) 


= k 1 (x + y 


Here 


(3.26) 


X .— Tp 0 t L gge — bfl T hfe L 4 L z L m ; x , 

is the combination appearing in the stress tensor, T m m = — 2X , and we define the useful 
quantity 

y ■= (g ~ 1 )Vb + (1 + p)L z + (1 + r)L m ix . (3.27) 

Notice that y involves only terms from the vortex lagrangian and vanishes identically in 
the scale-invariant case, for which p = r = — 1 and q = 1, while X is most usefully split 
between bulk and localized contributions through X = X B + X\ oc where X B = V B — L A while 
Aloe := Vb — L z - It is similarly useful to write Z = Z B + Z\ oc with Z B defined as above and 

-^loc = 2^ _ •^'gm- 

Once restricted to the symmetric configurations of interest the scalar equations simplify 
to 




and 


BW d 

Einstein equations 

The stress-energy tensor of the entire system now is 

1 



2 -v 2 


(3.28) 

(3.29) 


T M n = -w d M (j) d N <p + d M d N ip + e 2 ip 2 Z M Z N + e ^ A MP A N p + e p</, Z MP Z N p 

K Z 

+ 2 e r ^A M pZ N p + Z mp A n p ^J — g M N^L kin + ^gm + L pot + L gge ^ , (3.30) 
and so the nontrivial components of the matter stress-energy are given by (2.14), with X = 


Lp 0 t — Lgge (as in eq. (3.26)) while 

Q • — -I-'kin T -I-'gm T Lp 0 t + L 


gg e 


and 


Z ■— Lk; n L 


gm ■ 


(3.31) 


The lone nontrivial component of the trace-reversed Einstein equations governing the 
d-dimensional on-brane geometry therefore becomes 


R 


fc(d) ■- g^R-fiv - w2 + 


d 


W 2 BW d 


BW'W 


d -1 


= -2 kTX. 


(3.32) 
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The two components dictating the 2-dimensional transverse geometry similarly can be taken 
to be 


n {2) := g mn n. n 


R + d 


(W" B'W'\ 
\lT + BW ) 


= -K 2 X m m = -2 k 2 




(3.33) 


and 


B fW'\ 
W \~B~) 



(3.34) 


For later purposes a useful combination of these equations gives the ZD-dimensional Ricci 
scalar, 77m) = 77(<u + 77(2)- Given that the total lagrangian is given by L = L B + L v = 
(2 k 2 )~ 1 1Z(d' ) + g, where g is the total energy density, it turns out this can be written 


r _ _ 2X 

L - L b + Ly — — 

d 


(3.35) 


Yet another useful combination of the above equations is the (pp) component of the 
Einstein equations Q p p = —k 2 T p p which reads 

'W'\ 2 


2d 


B'W' 

BW 


+ ^ 2 +d ^ d ~ l ) 


w ) 


= 2k 2 (Z - X) . 


(3.36) 


This expression contains only first derivatives of metric and matter fields, and acts as a 
constraint on the solution as it is integrated along the proper distance p. 

Finally, we see from the above that the gauge fields only enter the Einstein equations 
through the combination L gge = L A + L z + L m j x = L A + L z , and so the Einstein equations 
are indifferent (unlike the dilaton field equation) to whether they are expressed using A MN 
or A mn . 


Control of approximations 

Since solutions to the classical field equations take up much of what follows, we first briefly 
digress to summarize the domain of validity of these solutions. The fundamental parameters 
of the problem are the gravitational constant, k; the coefficient of the bulk scalar potential, 
Vo (or Vo = 2g\l kA for 6D supergravity); the gauge couplings, e 2 (ip) = e 2 /A(ip) and g A (<p) = 
g A e^ 2 ; the scalar self-coupling, A(<^>) = Xe qv> , and the scalar vev v. To these must be added 
the dimensionless parameter, e, that measures the mixing strength for the two gauge fields. 
When discussing 6D supergravity we typically assume g A = g R and so can use g A and g R 
interchangeably. We also largely keep to the vortex parameter range A ~ e 2 . 

The energy density of the vortex turns out below to be of order e 2 u 4 and when A ~ e 2 
the transverse vortex proper radius is of order r v with 

fy = T 1 — • (3.37) 

ev 

The effective energy-per-unit-area of the vortex is therefore of order e 2 u 4 f 2 = v 2 . These 
energies give rise to ZD-dimensional curvatures within the vortex of order l/l 2 = K 2 e 2 v 4 and 
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integrated dimensional gravitational effects (like conical defect angles) of order k 2 v 2 . We 
work in a regime where nv <C 1 to ensure that the gravitational response to the energy 
density of the vortex is weak, and so defect angles are small and l v 3> r v . We also define the 
^-independent quantity r\, = 1/ev. 

By contrast, we have seen that far from the vortex the curvature scale in the bulk turns 
out to be proportional to l/£ 2 where 


£ = r B e 


K 

2 £«(</?) ' 


(3.38) 


Since our interest is in the regime where the vortex is much smaller than the transverse 
dimensions we throughout assume r v /£ <C 1 and so the parameter range of interest is 


< kv < 1. (3.39) 

e(¥>) 


As seen earlier, semiclassical reasoning also depends on the ambient value of the dilaton, 
<p, because it is e dlf ! 2 that counts loops in the bulk theory. Consequently we require 


e^<l 


(3.40) 


in order to work semiclassically within the bulk theory. But <p also governs the size of vortex 
couplings through A(y>) = \e qtf and e 2 (y?) = e 2 /A (cp) and we must check these remain small 
to trust semiclassical reasoning on the vortex. 


3.2 Dual formulation 

Because the gauge coupling to the vortex is magnetic, it can be useful to work with the 
Hodge dual of the Maxwell field A MN . In this section we restrict to the case of later interest, 
(. D,d ) = (6,4), though the same steps are easily also done for general d. 

The terms involving A M in the 6D action can be written 

L a + L mix + L lm = 1 e~^A MN A MN + | e r *Z MN A MN + i e MNPQRT V MNP d Q A RT , (3.41) 

where the functional integration over the 3-form lagrange multiplier, V M np , ensures A MN 
satisfies the Bianchi identity and so allows us to directly integrate A MN rather than the gauge 
potential, A M . Notice that because we wish the constraint dM( 2 ) = 0 also to hold on any 
boundaries we do not include a surface term to restrict the variation of V MNP there. 

The integration over A MN is gaussian and so can be performed directly, leaving V MNP as 
the dual field. Performing the gaussian integration requires an integration by parts, and so 
leaves a surface term 

Lst = +^V q (6 mnpqrt V mnp A rt ) , (3.42) 

to which we return later. The saddle point relates the 4-form field strength, = dV) 3 ), to 
the 2-form A( 2 ) as follows 

A mn = A mn + £ e (r+i)<A Z MN = - 1 e<P e MNPQRT F PQRT , (3.43) 

2*3. 
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which inverts to 


Fmnpq — +2 e ^ £mnpqrt-A r . 


The dual action (obtained after evaluating at this saddle point) then is 


(3.44) 


Tgge — Li z + L A + -^mix T L\ 

1 A (fh) Z, r ,rZ MN - - £ 


- L \(fk\ Z Z MN - " Jr+l)<j> MNPQRT z F 

— A Zj M N^ 0 /1 , £ t ^ MN r PQRT 


2-4! 


(3.45) 


+ 


^ p <t> i? 

2-4! PQRT 


F p QRT + L 


St 


L z + L blf + L F + L s t, 

where the last line dehnes L blf and L F , the latter of which also evaluates on shell to 

2 


L F : = 


1 


2-4! 


F, 


MNPQ 


= -i. = A(Q_ 

2 lVF 4 


pMNPQ = _L- 


o<t> 


(3.46) 


Keeping in mind that £ st = \f—g L s t does not change the bulk equations of motion, 
in these variables the contribution of the Maxwell held to the RHS of the dilaton equation, 
(3.25), becomes 

_ r)T _ . r A' .. 1 

(3.47) 


^2 d-^gge _ 


8(j> 


L z + (r + 1) L blf + L F 


instead of k 2 (—L A + pL z + rL m ; x ). Similarly, since C BLF = \f^g L blf is proportional to 
Z( 2 ) A F ( 4 ) it does not couple to the metric at all, so the stress energy coming from L gge is 


1 

3! 


rpMN _ ^ TzMABC TpN 

1 — — £ r ABC 


+ A Z MA Z N A ~ (L F + L z )g 


(3.48) 


and so contributes to X as A’gge = —^ T m m = —2 L z + (L F + L z ) = L F — L z . Furthermore, 
the 6D trace becomes T M M = 8 L F + 4L Z — 6 {L f + L z ) = 2 (L f — L z ) and so the full 6D 
traced Einstein equation is 

77 


(3.49) 


9 — + Lg m J 3L pot + L F L z 

— 2 Tkin T Tjgm T Tpot T L z L F ) I Lp 0t T L F L 


If we define y so that F\(f> = n 2 (X + y) remains true, then we are led to replace (3.27) 


with 


y — {q ~ 1 )Vb + ( 1 + ) L z + (r + 1) L l 


(3.50) 


As expected, L F drops out of this since the bulk Maxwell action does not break the scale 
symmetry. 

Trading A^) for F^\ in the surface term, L st , of (3.42) allows it to be written 

L* = +| V Q (e mnpqrt V mnp A rt ) = V Q (V MNP e* F MNPQ ) , (3.51) 
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where the second equality defines 


Fmnpq '■— F M npq 

Because the 4-form field equations imply V M 
solution gives 


- p r( t > f = 7, RT 

2 e t mnpqrt^ 

(i e <t>p MNp Q ) = o, evaluating L st 


L, 


St 


1 

' on—shell 4! 

and so on-shell the gauge action evaluates to 


= ~^e*F MNPQ F MNPQ = -2 L f - L 


BLF j 


L 


gge 


on—shell 


— L z + L blf + L F + L s t — L z — L F , 


(3.52) 
at a 4-form 

(3.53) 


(3.54) 


in agreement with the expected value, L gge = L z + L A , in the original variables. 

Although this dual formulation is equivalent to the original one, it makes several features 
usefully manifest. First, because F^ turns out to be proportional to *A( 2 ) rather than *A( 2 ), 
it provides a natural way to express the change of variables from A to A directly in the action 
rather than the field equations, even for nontrivial dilaton profiles. In so doing it generates 
the same renormalization of the Z kinetic term obtained earlier [4], Second, because the 
V — Z coupling term has the form of F^ A Zr 2 ) it is immediate that this term is independent 
of the metric and so does not directly gravitate. This can also be understood in the original 
variables [4] , in terms of a cancellation of localized contributions between L m [ x and L A . 

For the maximally symmetric configurations described above, evaluating F^ using the 
solution to the A( 2 ) held equation gives 

F, w \ k = e~^e^ u \ Kp gA pe = Q , (3.55) 


where e pv \ p = F^/—g is the 4D volume form built from the metric g pu . Notice that in the 
scale-invariant case (where r = — 1) these definitions imply F pu \ K —y s 2 F pu x K under the scaling 
symmetry. 


3.3 Vortex solutions 

This section describes an isolated vortex within a much larger ambient bulk geometry. Our 
goal is to establish that the presence of the dilaton couplings need not destroy the localized 
vortex solutions — with exponentially falling solutions beyond the vortex radius, r v — familiar 
from the dilaton-free case. We also wish to relate the properties of the vortex to the asymptotic 
behaviour of the bulk fields and their derivatives outside of (but near to) the vortex itself, 
with a view to using these in the discussion of bulk solutions given earlier. 


Nielsen-Olesen vortices 

The isolated Abelian-Higgs system contains vortex solutions where the local gauge symmetry 
is relatively unbroken in a core region [1, 2], and the fields approach their vacuum solutions 
outside of this region. Although we consider a more complicated gravitating vortex sector 
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that is coupled to a bulk scalar, as in S v , explicit numerical construction shows this only 
weakly perturbs the form of the localized vortex solutions in the parameter range of interest 
(as expected). 

We work in a unitary gauge for which T = is real (and for which the gauge fields Z M 
and A m are singular at the origin). We demand (as usual) the vortex fields approach their 
vacuum values away from the vortex, corresponding to i/j —> v and Z M —> 0 far from the 
vortex core. Inside the vortex we have scalar boundary conditions ^(O) = 0 and t//(0) = 0 
in addition to gauge boundary condition Zg( 0) = n v /e. This second boundary condition is 
chosen so that Z-flux quantization is satisfied within the vortex, 


Pv 


dZ = 2ir J d pd p Z e = 2ir Z e (p v ) - Z g ( 0) 

P=Pv o 


27m,, 


(3.56) 


where p v is a point chosen sufficiently far from the vortex core that we can assume the gauge 
field takes on the vacuum value Zq = 0, and the integer n v is the flux quantum of the vortex 
source in the region X v defined by 0 < p < p v . 

It is convenient when solving the vortex field equations to scale out the field dimensions 
by defining 

Zq = -1—— and i/j = vF(p). (3.57) 

e 

In terms of these the boundary conditions become .F(O) = 0 and F —>• 1 far from the vortex, 
while P(p) decreases from P( 0) = 1 at the vortex core to its asymptotic value P —> 0. In 
these variables, the vortex field equations read as follows 


1 

BW d 


( Bw d F 


nlP 2 F 
B 2 


+ A v 2 e q(t> F{F 2 


1 ), 


(3.58) 


and 


1 

BW d 




W d p' 

B 


+ (r + l)e (r+1) ^ 




e 2 v 2 F 2 P 
B 2 


(3.59) 


Examples of numerical solutions for the gravitating vortex profiles with dilaton interactions 
are shown in Fig. 3 and strongly resemble the nongravitating vortex solutions found in [4] with 
constant dilaton in the vortex region cj) = cf> v . In particular, they approach their asymptotic 
values exponentially over scales controlled by the ^-dependent masses m 2 z = e 2 ((j) v )v 2 and 
m\ = 2A (4> v )v 2 with 

® 2 ^) := A70) and A(0):=Ae 90 . (3. 

We wish to trace how the vortex solutions depend on cj), and to do so it is useful to rescale 
factors of A(^) into the coordinate p (or, equivalently, the transverse metric), by writing 


g mn dy m dy n = A (cj))r 2 (dp 2 + B 2 dd 2 ) = r 2 (</>) (dp 2 + B 2 d6 2 ) . (3.61) 
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With this choice the Ape equation becomes 


d, 


W d e~^Ap 0 


which has solution 


P BA{<f>) 

QBA(</>)e* 


= 0 , 


Ape — 


W d 


with integration constant Q = Qr 2 . The vortex field equations rewrite as 


BW d 9p 


BW d dpF 


(^) f (f 2 - 1 ) 


(3.62) 


(3.63) 


(3.64) 


and 


i (w d dpP\ 

BW d p \ B ) 


+ (r + l)e (r+1)< ^ 




F 2 P 


(3.65) 


Now comes the main point. To the extent that the gauge fields do not mix [ie e = 0) or 
that the dilaton is approximately constant in the vortex region (dp<j) — 0) these equations show 
that the vortex system depends on </> only through the one single, ^-dependent parameter 


'm ■■= (^) = ^ 

Furthermore, if p = 2r + 1 (as required if the Z M kinetic energy is bounded below for all <j >) 
then the ^-dependence of A(< j>) = e p ^(l— e 2 ) is simple and it is possible to make /3 independent 
of (f) by choosing p + q = 0. 

How important for these statements is the assumption that (f> not vary across the vortex? 
Depedence on <j) can enter if e is nonzero and 4> actually does vary across the vortex, as can 
be seen from eq. (3.59) or (3.65). But because (f>' ~ k 2 v 2 generically and eQ = er 2 Q ~ 
(e/ '9A){j‘v/r B ) 2 the vortex system is only weakly sensitive to the value of (j> for the one- 
parameter family p = —q = 2r + 1, since the ^-dependence arising from the cf>' term of (3.59) 
is of relative size 

dp<f>(r + l)ee r 2 Q e^ r+l ^ ~ K 2 v 2 (r + 1) (e^ r+1 ^ ~ (r + 1) ( e ( r+1 ^ , (3.67) 

\9AriJ \eg A J 

where the second estimate follows from the definition r v = l/ev and r B ~ n/g R . Because we 
are assuming 2r + 1 = p, the ^-dependence can be absorbed by the hatted couplings to give 

dp<j)(r + l)eer 2 Qe^ +1 ^ ~ (r + 1) (, (3.68) 

\eg A J 

which is suppressed in the parameter range g A /e ~ g R /e <C k 2 v 2 <C 1 required to control the 
semiclassical approximation and to ensure £ r v . 
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Figure 3: These plots demonstrate that the dilaton couplings do not ruin the existence of localized 
vortex solutions. The left plot contains a comparison of the vortex profiles F and P for a non¬ 
gravitating vortex in flat space (dashed curves) and a gravitating vortex coupled to the dilaton (solid 
curves). The right plot shows the dilaton profile in the vortex region for both cases. The gravitating 
dilaton is slowly varying in the vortex region, with the change of 0 over the vortex region being 
controlled by Ae 1 ^ ~ n 2 v 2 . The parameters used to generate the gravitating profiles are d = 4, 
£ = 0.3, p = 3, kv = 0.5, <(>(0) = 0, Q = 0.01 ev 2 , Vq = Q 2 /2 and the vortex sector is coupled to 
the dilaton through the generic choices (p, q, r) = (—1, 0, —1). The flat space profiles are generated by 
choosing kv = 0 instead. 


BPS special case 

In the special case where W = W v and 0 = </> v are constant in the vortex (for which a 
coordinate rescaling allows the choice W v = 1), then the vortex field equations are the same 
as apply in the absence of the dilaton once we make the replacement e —> e and A —> A with 
e 2 := e 2 /A(0 v ) and A := \ e q ^'. The vortex held equations in this case boil down to 


1 (P'Y _e 2 v 2 F 2 P 

b\1b)~ W~ 


while the 0 equation becomes 


1 

B 



n 2 P 2 F 

B 2 


+ \v 2 F (F 2 - 1) . 


(3.69) 


(3.70) 


The solutions to these equations are particularly simple when e 2 = 2A, since then 
eqs. (3.69) and (3.70) are equivalent to the first-order equations, 


BF' = n v FP and 


n v P' 

~Jb 



(3.71) 


We show later that W = 1 and 0 = 0 V also solve the bulk held equations when e 2 = 2A, and 
so this choice provides a consistent solution to all the held equations. Such solutions naturally 
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arise when the vortex sector is itself also supersymmetric, since supersymmetry can require 
e 2 = 2A and the vortices leave some supersymmetry unbroken. 

When eqs. (3.71) as well as (f> = <p v and W = 1 hold, they also imply 

L kin = ^ W) 2 = y VZ m Z m = L gm , (3.72) 

and 

L z := i K(fo)Z mn Z mn = \(V~ v 2 ) 2 = V b , (3.73) 

which further imply that the vortex-localized contributions to Z and X cancel out: Z\ oc = 0 
and X\ oc = 0. But the bulk contribution to Z also vanishes if <f>' = 0 and — as can be 
seen from eq. (3.34) - it is the vanishing of Z that allows constant W to solve the Einstein 
equations. Finally, the dilaton field equation with constant (j) requires X B +y = 0 everywhere, 
and so separately evaluating in the bulk and vortex implies X B = y = 0 separately. Although 
y = 0 can be ensured using the scale-invariant choices p = r = —1 and q = 1, vanishing X B 
in general either requires a condition on the bulk gauge held, Q (which need not agree with 
what is required by flux quantization) or a runaway to —> 0. 

Finally, the vortex part of the action evaluates in this case to the simple result 


S - J_ 


pv 


d 2 yV~g (hoc = 27T / d pB 

x v o 


L<h + V b + Lz 


= 7 m v v 


(3.74) 


where the second equality also defines the localized energy density £q oc = L\ oc . 


3.4 Integral relations 

In this section we generalize the integral relations described earlier for the bulk system to 
include the vortex sources. Instead of integrating only over the exterior region, B ex t, we now 
instead integrate over the small regions, X±, containing each vortex source, and thereby learn 
how the vortex determines the boundary conditions on the interface with B ex t. Using these 
boundary conditions for the integrated bulk solution is equivalent to integrating the bulk- 
vortex field equations over the entire space, X tot := B ex t U X + U A_, which is smooth and 
compact and without boundary. In what follows we generically represent by X v = {X+, A"_} 
when we do not need to specify which source is of interest. 


Integration over near-source pillboxes 

We first integrate just over X v to find the UV theory’s perspective on the general matching 
conditions [8, 9] that relate near-source derivatives to properties of the source. These are 
useful in developing the effective theory of the next section that treats the vortices as point¬ 
like sources, or branes. 
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Figure 4: These plots show the numerical bulk solution for n = 1 BPS vortex sources (ft = 1) with 
scale invariant couplings to the dilaton (p, q, r) = (—1,1, —1) and kv = 0.4. For these choices, the field 
equations are solved by a constant dilaton and warp factor: W = 1 and </> = 0. Because ft = W' = 0, 
this solution falls into the simple class of rugby ball solutions described in §2.3. In the top left plot 
the metric function B is plotted against the same solution for a sphere of proper radius £ = 2507rr v , 
B S phere = £ sin(p/£). The vortices (which cannot be resolved in these plots) introduce a defect angle 
into the bulk metric such that B = a£sin (p/t) in the bulk. The defect angle, a, can be determined 
by extrapolating B' to the apparent singular point at p* ~ 0. This yields 1 — a ~ 0.08 = k 2 T/2tt = 
(kv) 2 / 2 as expected from (3.74) and (4.25). 


Maxwell fields 


We start with the gauge-field equations. The simplest of these to solve is the Maxwell equa¬ 
tion, (3.14), since this does not depend directly on the fields Z M or if. The solution is as 
before 



QBe^ 
W d ’ 


(3.75) 


where Q is an integration constant. This enters into the Einstein equations, (3.32), (3.33) 
and (3.34), through the combination L A = \(Q/W d ) 2 e^. 

Suppose now we take a test charge that couples only to A M and ask how much flux it 
measures when taken around the vortex. We do so by moving around the edge of V, remaining 
everywhere outside the vortex. The edge of the pillbox is at a distance p v > r v from the vortex 
so vortex fields are exponentially small, but we also choose p v <^ £ so it contains a negligible 
fraction of the external bulk. The flux seen by this charge is 
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*a(*v) : 


= f A ( 2 ) = f ( A m - ^ r+l) * z w) 


= 2ir 


pv Pv 

Q J dp -e j dp e^ r+1)<t> Z pe 
. o o 


(3.76) 


The first term in the last equality gives the amount of bulk flux lying within p < p v in the 
absence of the vortex source, and so is negligibly small in the limit r v —> 0. The second 
term does survive this limit, however, because even though bulk fields vary slowly over the 
small vortex volume, the profile for Z p g is strongly peaked in such a way as to give the total 
quantized Z-flux, 

= (3.77) 

X v 

So eq. (3.76) shows that the mixing of A and Z gauge fields through L m ; x implies the test 
charge sees a vortex-localized component of flux despite it not coupling directly to the Z field 

$ A (X V ) ~ -e $ Z (X V ) e (r+1)< ^ ~ 2 ^- e (r+1)< ^ , (3.78) 

e 

where the approximation is true to the extent that f v <C l. This localizes part of the external 
A flux onto the vortex. 


Dilaton 

Integrating the dilaton field equation, (3.29), over a vortex-containing pillbox gives 

[ bwW) p=pv = ^f d 2 y ^ w d (x + y ) = ^x + y ) v , (3.79) 

x v 


where we use that BW d (j)' vanishes at the vortex centre p = 0. We also repurpose the angle 
bracket notation of (2.24) to more general integration regions, X v . This exact result expresses 
how the near-source limit of q 7 just outside the vortex is determined by the detailed vortex 
profiles (keeping in mind that for both regions X±, (f)' in this expression is evaluated for a 
proper-distance coordinate for which p increases as one moves away from the vortex). 

If p v lies within the Kasner regime — for which (2.30) applies — the the left-hand side 
of (3.79) becomes 



^ zB 0 W d ^ 


\ dw-\-b— 1 


l + O 



zT 


l + O 



(3.80) 


where p v := p v — p* and the second equality defines the quantity T := BoW d /£ and uses 
the linear Kasner relation dw + 6 = 1. When combined with (3.79) this shows how vortex 
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properties constrain combinations of bulk parameters (such as z) not already fixed by the 
bulk held equations. 

Metric 


Similar conditions are obtained by integrating the Einstein equations over X v . The trace- 
reversed Einstein equation, (3.32), governing the curvature 7 Z^) integrates to give 


B (w d ^j 

which uses the boundary condition B (W d Y = 0 at p = 0 and rewrites [B (W d \\ Pv using 
(2.30). This relates a different combination of bulk parameters to vortex properties. Inte¬ 
grating the (99) trace-reversed equation instead implies 


= -^l A ( w ~ 2 ),+ 2k 2 <*>,]. (3.81) 


dwT 


1 + 0 l ^ 


bT 


1 + 0 







(3.82) 


because of the boundary condition B'W d = 1 at p = 0. 

Notice that the powers z, w and b are not independent since the bulk held equations 
imply they must satisfy the Kasner conditions (2.31), and so the right-hand sides of the 
above expressions also cannot be completely independent in the limit p v <C l. The resulting 
relations among the vortex integrals are developed in more detail in §4 and play an important 
role in determining the off-brane components of the bulk stress energy in the effective theory 
applying at scales where the vortex size cannot be resolved. Because these relations follow 
from the bulk Einstein equations they can be regarded as general consequences of stress-energy 
conservation for the vortex integrals. 


Integration over the entire transverse space 

We see that the integral relations of §2 give bulk properties in terms of near-vortex deriva¬ 
tives of bulk helds, and the integral relations just described then relate these near-vortex 
derivatives to explicit vortex integrals. The resulting relation between bulk properties and 
vortex integrals is more directly obtained by integrating over all of the transverse dimensions, 
Wot = B ex t U W|_ U X-, at once. In such an integral all boundary terms cancel, as they must 
for any smooth compact transverse space. In the case of the Maxwell field integrating over 
the entire transverse space gives the flux-quantization condition, as discussed earlier. 

Dilaton 

Integrating the dilaton field equation, (3.29), over the entire compact transverse dimension 
gives 

(x + y) M = 0 « (X.) M + E(-*ioc + y\ , (3.83) 

V 

where the approximate equality drops exponentially suppressed vortex terms when replacing 
a localized integral over the entire space with a localized integral over the source regions. 
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Integration over the transverse space can be regarded as projecting the held equations onto 
the zero mode in these directions, and so (3.83) can be interpreted as the equation that 
determines the value of the dilaton zero-mode. (This conclusion is also shown more explicitly 
from the point of view of the effective d-dimensional theory in a forthcoming analysis [33].) In 
the absence of the sources this zero mode is an exact hat direction of the classical equations 
associated with the scale invariance of the bulk held equations (for instance X B = 0 for the 
source-free Salam-Sezgin solution [34]) and the vortex contribution to (3.83) expresses how 
this hat direction becomes hxed when the sources are not scale-invariant. 

Metric 

Integrating the trace-reversed Einstein equation over the entire transverse space leads to 

H-HHoT 0 ' 

and 

« = - 2 « 2 <r> tot = d^(L ) ta , (3.85) 

which shows how it is the stress-energy transverse to the source, (T) to t = — |(T)tot) that 
ultimately controls the size of the on-source curvature [7]. Using (3.83) to eliminate X from 
the right-hand-side similarly gives 

^A,„, = 2 « 2 W«“I> 2 Mv. 

V 

whose approximation drops exponentially suppressed terms. 

A hnal useful rewriting of these expressions uses the formula relating the D- and d- 
dimensional gravitational couplings, k 2 and k 2 d respectively. Dimensionally reducing the 
Einstein-Hilbert action shows that this states 

i = 4 A -2 >.„, • < 3 - 87 ) 

K d K 

when the would-be zero-mode for 4> is evaluated at the solution to its field equations, and so 
eqs. (3.84), (3.85) and (3.86) at this point become 

R = >#)„ = -2*1 <*}„ = - (^ 2 ) 4 <e - ■ (3-88) 


(3.86) 


(3.84) 


4. Sources - effective IR description 

This section takes the point of view of a low-energy observer, and recasts the expressions for 
the UV theory found above into the language of the effective field theory appropriate in D 
dimensions at scales much larger than the transverse vortex size, r v , but smaller than or of 
order the KK scale, l. We specialize for concreteness’ sake to the case (D, d) = (6,4), though 
our conclusions also hold more for general D = d + 2. 
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4.1 The EFT with point sources 

If the relevant length scale of an observable r Q b s exceeds the length scale of the vortex, 
fobs f v , then effects of the vortex can be organized as a series in the small quantity r v /r Q b s . 
For sufficiently large r 0 b s , the internal structure of the vortex becomes irrelevant and it can 
be replaced with an idealized point-like object. This is why Abelian-Higgs vortices without 
localized flux are well-described by the Nambu-Goto string action at long distances [1, 2 ], 
We here generalize this to include brane-localized flux, extending [4] to include dilaton 
dependence. It is most convenient to do so using the dual formulation of the bulk action, 


S B = - d 6 xy/^g 


2 g 


^2 g MN (TZ-mn + d M (j)d N (j)) + -— (Fmnpq) + e ' 


(4.1) 


and to include localized flux in the brane action we include the first subdominant term in a 
derivative expansion 10 

S e s = ~ ^ [ d 4 c \/—7 

v Jx=Zv(<t) 

(4.2) 

where 7 ^ 1 /(c) = gMNd^z^ 1 d v z™ is the induced metric at the position of the brane (with z™(cr) 
denoting the brane position fields) and £^ uX p is the totally antisymmetric 4-tensor associated 
with this metric. Since in what follows our interest is not in the dynamics of the brane 
position modes we assume a static vortex and choose coordinates so that it is located at fixed 
y m = VT and identify a 1 * = x^ so 7 ^( 2 ;) = g^, v (x,y v ) = W 2 (y,)g^{x). It is clear in both 
the UV and IR theories that the term linear in F (4) does not gravitate because it is metric 
independent, though this can also be inferred using the dual variables as in the dilaton-free 
case [4], 

Because the effective theory cannot resolve the vortex structure it also cannot distinguish 
between the quantities p v and associated p* used in previous sections. For each vortex we 
define the brane position in the effective theory to be the corresponding place where the 
external metric is singular when extrapolated using only bulk field equations. In practice this 
situates them at p = p* in the coordinates used earlier, where p* = {p+, p~} is determined by 
B{p ic ) = 0 , so we use the notations </>* = </>(p*) = 4> v = 4>(y v ) interchangeably for the various 
bulk fields. 

We show in this section that captures all of the physics of the full vortex action, 
up to linear order in the hierarchy r v /t, provided that the parameters T v ((j>) and £„((/>) are 
chosen appropriately. In general, these can be ^-dependent quantities and we identify this 
dependence by demanding agreement between the predictions of and S v to this order. 
Working to quadratic or higher order in f\,/£ would require also including higher-derivative 
terms in To connect the point-brane action to the bulk fields we promote it to higher 

10 The check here conforms with the notation of [4] and is to remind that the use of the 4-form field auto¬ 
matically unmixes the gauge kinetic terms. 




-=E 


d 4 x £ v =: / d 4 x C eS , 
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dimensions using a ‘localization’ delta-function, S(y). More precisely, we write 


s* = Ej i °* £ ’( s ^r)’ < 4 - 3 > 

where 5(y) is localized around zero and it is normalized so that integrating over a single source 
region X v gives J x d 2 y6(y — y v ) = 1. Performing the integration over the extra dimensions 
recovers the brane action in (4.2). 

As in [4] a key question asks what fields the localizing function depends on since this 
affects how 5 e ff enters the bulk field equations. Although we assume in what follows that 
8(y) is independent of the fields A M , <f> and g pu , we cannot also do so for the metric components 
9mn because it is designed to discriminate points based on proper distance from the vortex 
center. But because this metric dependence is only implicit it complicates the calculation 
of the brane’s stress-energy components Tmn ■ One of the purposes of this section is to show 
how to determine these components without making ad-hoc assumptions about how 8(y) 
depends on g mn . instead deducing them using properties of the bulk Einstein equations. Our 
conclusion ultimately is that the naive treatment of ignoring metric dependence in 5(y) need 
not be justified in the presence brane-dilaton couplings, and a cleaner way of inferring how 
a brane interfaces with the bulk is provided by directly relating the near-brane boundary 
conditions of bulk fields to derivatives of £ e ff (along the lines of [8, 9]). 

4.2 Parameter matching 

We start by matching the coefficients 11 T and £ by comparing with the UV theory. A 
direct way to do so is by dimensionally reducing the UV action, and we verify that this 
also is what is required for 5’ e ff to agree with S v for observables like the components of the 
stress energy. When doing so it is crucial to notice that any such a comparison between 
the UV and IR theories need only be done up to linear order in f v /t since it is only at 
this accuracy that the action (4.2) must capture the physics of earlier sections. This allows 
considerable simplification since integration of any slowly varying bulk quantity over X v 
vanishes quadratically with r v for r v <C l, allowing any such terms to be dropped to the 
accuracy with which we work. In what follows we accordingly take the formal limit r v —> 
0 when discussing such integrals, by which we mean we drop terms that vanish at least 
quadratically in f v /£ in this limit. We take care not to similarly drop terms suppressed only 
by a single power of r v /£, however. 

The term linear in F pv \ p in the UV action is given in (3.45) as 

L blf = ~^e( r+1 ^e MNPQRT Z MN F PQRT = -^ e^ +1 ^e mn ^Z mn F puXp , (4.4) 

and because the 4-form Bianchi identity, dF( 4 ) = 0, ensures the components F pu \ p cannot 
depend on the transverse coordinates, y m , we know F pu \ p cannot be strongly peaked (unlike 

^Although we often drop for simplicity the subscript V from T, £ and p* the reader should keep in mind 
that these quantities all can differ for different branes. 
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Z mn ) iii the off-brane directions. Consequently integrating over the vortex area and comparing 
with the corresponding term in the IR theory gives 

CvOM = J lim [ d 2 y B e ( r+1 >* e mn Z mn ~ - ( e (r+i)<t>. (4.5) 

2 f v ->o J \ e J 

x v 

where the approximate equality neglects the small variations of <fr away from a constant value, 
(j) v , at the vortex position, and uses Z-fiux quantization to evaluate $> Z (X V ) = 3 / d 2 y e mn Z mn = 
—2irn v /e. 

To match the tension T we compute the source contribution to the energy density. In 
the IR theory we have 


Tfiv — 9/if 


V B + Ld 



(4.6) 


and this is to be compared with the localized contribution to the stress energy of the UV 
theory integrated across the vortex, 

T llv = i F^ m - n F u m .. n - g, iv (V B + L ( j, + L F + L z + L^ + V b ) 

= ~9ta> {Vb + — L F + L z + L^i + 

= -giif(j>B + eioc) ■ ( 4 -7) 


Comparing T^u with lim^ i ,^o(T /liy ) v reveals the localized contribution to the energy density in 
the UV theory to be (temporarily returning to general d) 


W* 4 T[<Kp*)] = lim 


d 2 yBW 4 


L z + L\fy + V b 


x v 


lim (Lz + Ly + Vb) = lim (g\ oc )v ■ (4.8) 

r v —>0 \ iv f r —>0 


We pause here to highlight one important feature of this result. The near-brane behaviour 
of the warp factor is W oc p w and for w > 0 the quantity IT 4 = VU 4 (p*) formally vanishes. 
The power law vanishing of W can be reinterpreted as the logarithmic divergence of the 
field c = ln(W) and such divergences are common in theories with higher codimension brane 
sources. These divergences can be classically renormalized into the brane couplings [10, 15] 
and in this case the tension would be renormalized such that the physical combination W^T 
remains finite. The UV complete theory provides an explicit regularization of this divergence, 
since the vortex physics intervenes near the source to ensure W > 0 everywhere in the vortex 
region. 

Unlike for £, the result in (4.8) gives the (^-dependence of T only implicitly, so we next 
display this dependence more explicitly. We first compute how T depends on cj) assuming 4> 
to be constant over a vortex. We expect the errors we make by doing so to be suppressed 
by powers of r v /£, and come back to verify this estimate shortly. The ^-dependence of the 
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tension is determined by the ^-dependence of vortex integrals like 


{Vb) v = \J d2 y BW 4 (V ’ 2 - A 2 , 

x v 

(L gm ) v = \J d2 y BW 4 e\^ 2 Z m Z m , 
x v 


(L z ) v = i J d 2 y BW 4 A($)Z mn Z mn 
x v 

(l* kin) v = \J d2 y Bw A d m ^d m ^, (4.9) 

Xv 


and earlier sections show that the ^-dependence of the vortex prohles appearing in these 
are fairly simple once expressed in terms of dimensionless variables, F and P. Then the 
implicit ^-dependence within the profiles themselves arises only through the combination 
j3((j)) = 2\(c/))/e 2 (4>) , with A ((/>) := \e q ^ and e 2 /e 2 (c j>) := A(<f>) = e p ^ — e 2 e^ 2r+1 ^. We now 
ask whether any additional (^-dependence arises from the integrations to set the scale of the 
above integrals. 

To this end, it is useful to return to the variables used in (3.61) in which B is dimension¬ 
less, as is the radial coordinate p. The integration measure d 2 y BW 4 = (1/eu) 2 d 2 y BW 4 and 
we have, for example, 


/d 2 sw 4 (f 2 -1 ) 2 — d J^yBw* (F 2 -iy 

x v x v 


Similarly 



(4.10) 


(4.11) 


while 

{Lp ki„) v ^y J d 2 y BW 4 (F') 2 , (4.12) 

Xv 

and so on. In these expressions the integrands are all proportional to e 2 v 4 , but the e 2 
dependence cancels the factors of r 2 = l/(e 2 v 2 ) coming from the integration measure to leave 
integrated results that again depend on 0 only through /3(0) oc e q ^A(cj)). 

Consequently 

T(</>) = t0((/>)] , (4.13) 

and in particular T is ^-independent for the one-parameter family of choices p = —q = 2r + l. 
Because C(^) cx e^ r+1 ^ having both T and £ be ^-independent happens only in the special 
case of scale invariance, for which p = r = — q = 1. This inference of the ^-independence of 
the tension is verified numerically, as seen in Fig. 5. 

What happens once we drop the assumption that <p' is negligible within the vortex? In 
this case our earlier estimate — eq. (3.68) — of the leading (^-dependence of vortex prohles 
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p/f'v pl r v 

Figure 5: A demonstration of the ^-independence of T v . Two solutions are presented for the de¬ 
coupling choice p = —q = 2r + 1 = 2 that differ only in the value of the dilaton in the vortex 
<f> v ~ ^(0) = {0,-2}. The light (grey) lines represent the solution for the choice <)>(0) = —2 and the 
blue lines represent the solution for </>( 0) = 0. While the physical mass scales that control the size 
of the profiles, m%(4> v ) = 2Au 2 e ? A and m|(</> v ) = e~ p ^ v e 2 v 2 /(1 — e 2 ), are demonstrably heavier for 
(j> v ~ —2 (since these profiles fall off much faster), the defect angle [TT d i?'] ~ 1 — that measures 
the tension of the vortex is independent of <f> v . The other parameters of this solution are d = 4, e = 0.6, 
/3 = 0.8, kv = 0.4, Q = 1.6 x 10 " 4 and V 0 = Q 2 / 2. 


implies that the leading ^-dependent corrections to quantities 
KK-suppressed form 

6T V _ 5(X loc ) v 

T v (diloc)v 

which uses g A ~ g R and p = 2r + 1 as well as r v ~ 1/ (ev) and i 


like T v and (Ai oc ) v 

(r + 1 )kv , 

~ n/g R - 


have the 

(4.14) 


4.3 Near-source matching conditions 

As noted earlier the hidden dependence of the localization function 5(y) on g mn complicates 
the inference of how branes contribute to the bulk field equations. In this section we give a 
^-independent way of expressing the bulk-brane interaction wherein knowledge of the brane 
action, S e s, directly gives the near-brane asymptotic derivatives of the bulk fields. This 
connection is the analogue in the IR theory of the expressions in §3.4 relating the near-vortex 
derivatives of bulk fields to integrals over the vortex. The agreement between the UV and IR 
descriptions of these boundary conditions provides a check on the matching between the two 
versions of the theory. The discussion found here also parallels the dilaton-free case [4] and 
this source can be consulted for more details. 

Dilaton 

Since the general discussion is simplest for a scalar field we treat the dilaton first. In the 


- 36 - 











effective theory of point-like branes the field equation for the dilaton reads 


1 


n<f> = v B 





(4.15) 


As before we integrate this equation over a source-containing pillbox to isolate the near-brane 
derivative of the dilaton, recognizing that this pillbox can be taken to have infinitesimal size 
in the effective theory (within which the vortex size cannot be resolved). On the left-hand 
side the result is 

lira (D0) v = 2 t r lim ( BW ( V') = 2vr zT, (4.16) 

where the last equality uses the near-brane asymptotic bulk solution (2.30) and again uses 
the definition T := BqWq /£. Performing the same operation on the right hand side gives 


lim(v B -L a + = W\p*) (w) - _ , (4.17) 

where the smoothness of the bulk sources ensures their integral does not survive the limit 
r v —> 0. The localized sources do survive this limit, however, and give the final result 


zT = 



k 2 W° 

2 tt 


t“M) - Vi= 

' P* 


^ (SSes\ 

2vrV^ V S( P ) P=Pi , 


(4.18) 

in agreement with [8, 9]. Eq. (4.18) is useful because it directly extracts the impact of the 
brane-dilaton coupling on the bulk dilaton without reference to the localization function 5(y). 

This argument shows it is useful to divide quantities like g, X and Z into a localized 
piece, whose integral survives the point-like limit, and a smooth bulk piece that does not. 
For instance: X = X B + X\ oc , where the bulk part, X B := V B — La , depends only on bulk 
fields and so satisfies ( X B ) V —> 0 as f v —> 0. The same need not be true for vortex-localized 
quantities like X\ oc = 14 — L z , y and so on. 

The analogous relation between the near-vortex derivative of the dilaton and the vortex 
sources in the UV theory is given by (3.79), which we rewrite for convenience of comparison 
here 


lim BW d $ 
p->p* 


im (Aioc + y ) = 






lim 


A' - 

nVu -I- T,r 


(4.19) 


Comparing with this UV version allows a check the consistency of our inference of the (f>- 
dependence of L e ff. Comparing terms linear in gives 




(r + l) lim - / d 2 yBe^ r+1 ^e mn Z n 
f v —*o 2 I 


x r 


(4.20) 


which is consistent with the earlier result (4.5), assuming the dilaton is approximately constant 
in the region X v and provided differentiation with respect to cj) is performed with fixed vortex 
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fields. Similarly comparing the 4-form-independent terms gives 

(w d f;) =lim (qV b +^L z ) , (4.21) 

V / p=p* r v ->0 \ A / v 

which is also consistent with the earlier expression (4.8). 

Metric 

A similar argument relates near-source metric derivatives to properties of the brane action. 
As above one integrates the Einstein equations over a region X v enclosing the vortex and 
finds two kinds of terms that survive the limit r v —> 0 of vanishingly small vortex size. One 
such class of terms comes from the vortex parts of the stress energy while the other come 
from terms inside the Einstein tensor involving second derivatives with respect to p, with all 
other contributions not singular enough to survive the small-vortex limit. 

The simplest equation to analyze is the (pp) Einstein equation, (3.36), since this involves 
no second derivatives at all. Consequently its integral over X v simply states 


T% = lim (TP p ) v = 

r v —>■ 0 


lim {Z\ oc 

TWO 


Aloc)v ^ 0 . 


(4.22) 


Physically, this is a consequence of dynamical equilibrium for the non-gravitational micro¬ 
physics of which the vortex is built, since this requires there to be no net radial pressure. 

For all of the remaining Einstein equations the Einstein tensor does include second deriva¬ 
tives and so their integration over X v leads to the following relation between the near-source 
derivatives of the metric and the metric derivative of the source action 


27t lim 


> 0 L 


sJXTg (jpj _ KgV)] Pv = -n 2 ^jxr g pi = -2 k- 


J o 


,(55eff 

39ij 


(4.23) 


where the derivative on the right-hand side is with respect to the metric evaluated at the brane 
position and is the extrinsic curvature for surfaces of constant p (with K = g 13 Kij). The 
indices i and j run over all coordinates but p. 

The derivative on the right-hand-side can be taken reliably for on-brane components of 
the metric (ij) = {pv) using the action (4.2) and gives 


QQiiv * 


For the metric ansatz of interest the extrinsic curvature components are K^ u = WW'g^ and 
K ee = BB ', so the on-brane components of (4.23) give 


lim 

r v —>0 


1 - BW 4 


, , ,, W' B' 

{d ~ 1] W + ~B 


P=Pv 


lim 1 

p->p* 


d- 1 


B 



B'W d 


k 2 W 4 T 

2tt 

(4.25) 
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This result is to be compared with the appropriate linear combination of (3.82) and (3.81) in 
the UV theory, keeping only those vortex-localized terms that survive in the limit f v —>• 0: 


lim 

p—»p* 



d- 1 
d 


B 



B'W° 


K 2 

- x- lim {q\ oc + X\ oc 


K, 


(4.26) 


The last equality here uses (4.22), leaving a result consistent with our earlier identification in 
eq. (4.8) that W d T v ~ lim^-xj^io^v 


4.4 Brane and vortex constraints 


We now turn to the remaining case, where (ij ) = (99) in (4.23), which gives 


lim 

r v —>0 


B 



Pv 


lim B ( W d 
p—tp* V 




9ee 



(4.27) 


The trouble with this expression is that the dependence of S e s on gee is only known implicitly, 
so we cannot perform the differentiation on the right-hand side to learn about the near-brane 
derivatives on the left-hand side. Fortunately we may instead read this equation in the other 
direction: it tells us the right-hand side because the radial constraint - eq. (3.36) - already 
determines the derivatives on the left-hand side in terms of known quantities. It is ultimately 
this observation that allows us to determine (X\ oc ) v and (Z\ oc ) v separately in terms of the 
quantities T and ( [8, 9]. 

To see this we first compare with the corresponding UV expression, using (3.81) to rewrite 

lim B (W d )' ~ - — lim {X loc ) v , (4.28) 

p->p* V / 7r f v —»o v 


and so 


n~6 

I e 



lim(^° c + Z loc ) ^, 


(4.29) 


which confirms that the vortex stress-energy component T e e captures the integral of the 
vortex-localized part of T e q in the UV theory. 

To fix T e o we follow [8, 9] and again use the constraint equation (3.36), but rather than 
integrating it over the vortex we instead evaluate it at p = p v . just outside the vortex, and 
solve it for B (W d y to find 


B (W d ) = - ( B'W d ) + (B'W d ^ J 1 + E , (4.30) 

where we define 

/ JD \ 2 

E := f — j [2k 2 {Z b - X B ) - W- 2 R] . (4.31) 

Because this gives B ( W d y in terms of B'W d and other known quantities we use it in (4.27) 
to get an explicit expression for T e q (and so also for (Z\ oc ) v ~ (T) oc ) v ). 
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So far eq. (4.30) assumes only that p v is sufficiently far from the vortex that localized 
contributions to X and Z are exponentially small. However, since p v ~ r v in size we can also 
drop terms quadratically small in p v when comparing with the point-brane theory. Since both 
the R and X B terms in (4.31) are proportional to B 2 (p v ) ~ p 2 ~ f 2 they can be dropped in this 
limit. By contrast, the term containing 2 n 2 Z B {p v ) = [<j>'(p v )] 2 need not vanish quadratically 
in this limit, since the asymptotic power-law form (2.30) allowed in this region implies 


lim E ~ lim 2 k 2 ( — ] Z B ~ lim 
f v — s-o p-»p* \B ) p-*p* 




(4.32) 


and the final equality uses (3.80) and (3.82) to rewrite the right-hand side in terms of Kasner 
powers. Indeed, using this final result in (4.30) and trading the remaining terms for Kasner 
powers leads to 

(d - l)w ~ -b + 6^1 + ) (l) > ( 4 - 33 ) 

which reveals it not to be independent in this limit of the two Kasner conditions, (2.31). 

In passing we pause to remark on a point already alluded to in earlier sections: that the 
constraints imposed by the Einstein equations also imply the existence of relations among 
vortex integrals — like {X\ oc ) v and (T)v — f° r arbitrary vortex microphysics in the point- 
source limit. The constraint is found by eliminating B' and W' in terms of vortex integrals 
using the near-vortex expressions (3.82) and (3.81), leading for instance to 


lim E — 

r v —>0 



K 2 (*loc+y) v 


2t r(^l° c 2 Ti oc ( d d 1 ) ) V/ 


(4.34) 


In limit of small K 2 (y) v and n 2 (X\ oc ) v the resulting constraint simplifies to 


k 2 (Ti oc ) v ~ 


^(y)l 

47T 


(4.35) 


This brane constraint can be verified using explicit numerical solutions, as is illustrated in 
Fig. 6. The suppression it implies for k 2 {X\ oc ) v relative to n 2 {y) v is the analog of the suppres¬ 
sion found earlier, (2.33), between the various Kasner exponents that imply w and 1 — 6 are 
order z 2 when z«l. This is seen most explicitly by dividing (3.80) by (3.81) and evaluating 
in the r v —> 0 limit, which gives 


2 (Ti oc ) v ^ dw 
(dfloc +T) V ““ 


z 

2 


+ 0(z 2 ) 


(4.36) 


where the approximation follows from (2.33). The integrated vortex sources (X\ oc ) v and (y) v 
cannot adjust independently to leading order in 6 and r v /£. This is explored in more detail 
in the next subsection. 


- 40 - 









Figure 6: Two examples of the vortex constraint (4.35). The light (grey) line is calculated from the 
expression k 2 (T| oc ^ ~ — tt [B(W d )'] p and the dark (blue) line is calculated using using the second 
brane constraint. Once evaluated outside of the vortex p v > r v « r v , the two independently calculated 
quantities are in perfect agreement, and their failure to agree inside the region p v <, r v reflects the fact 
that the vortex constraint does not hold locally in the source. In the first plot the parameters are 
d = 4, e = 0, /3 = 3, kv = 0.5, 0 V ~ 0(0) = 0, Q = 2 x 10” 4 and V 0 = Q 2 / 2. The vortex is coupled to 
the dilaton via the choice (p, q) = (0, 2) and r is not relevant because gauge kinetic mixing has been 
shut off. In the second plot, the only change in parameters is that 0(0) = 1. 


To summarize, it is possible to be very explicit about the dependence on S e g of k?(X\ oc ) v 
and k 2 (Z\ oc ) v in the limit when these vortex integrals are small. In this case ( y) v is given by 
differentiating S e s by eqs. (4.18) and (4.19), 


zT = 



2 



~ — lim (T) , 

27Ttwo X /v ’ 


(4.37) 


Then ( Z\ oc ) v and (Ti oc ) v are obtained from (4.22) and (4.35). 

This section also shows how errors can arise when treating the localizing function 5{y ) 
as a naive delta-function that is independent of of g mn . Such a treatment would mistakenly 
conclude that the brane action in (4.2) is independent of the bulk metric, and give 


lim B(W d y = —-=—gee ( r eff ^ = 0 (naive result). (4.38) 

p^p* ny/-g \ Ogee J 

This not consistent with (4.28) unless the vortex source has vanishing (dq oc ) v in the limit 
r v —>• 0, which is not necessarily true, since the vortex must satisfy the simplified vortex 
constraint (4.35). In the next section we estimate this quantity’s (often nonvanishing) size. 


5. The scale of the response 

This section summarizes the implications of the previous sections for the generic size of back- 
reaction effects on physical quantities. In particular we broadly scope out how the size of 
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the on-vortex curvature, R, varies with the parameters p, q and r governing the size of the 
vortex-dilaton couplings. 

As (3.88) suggests, the d-dimensional curvature is sensitive to p, q and r through the size 
of K d(y) tot- We now trace how the size of R can drastically change as these parameters take 
various limits. The size of the quantity {X\ oc )v also varies strongly in these limits since it is 
also related to (T)tot ~ because of (4.22). 


5.1 Generic case 


Before examining special cases we first establish a baseline by considering the generic case. 
Recall from the definition of y in the dual variables, 


y — (q — 1)14 + ( 1 + -^- ) Lz + (j‘ + 1 )L blf , 


(5.1) 


that each term contains localized vortex fields and for generic choices of parameters these all 
integrate over a vortex to give 


CV)v ~ (f?ioc)v ~ T ~ zr , (5.2) 

much as we saw earlier. 

When this is true, it follows from the vortex constraints (4.35) and (4.22) that the vortex 
integrals of the transverse stress-energy components are suppressed relative to (y) by an 
additional factor of k 2 v 2 , 

k 2 (Aoc) v - k 2 (A’ 1oc ) v ~ kV . (5.3) 

This suppression is perhaps not surprising given that these quantities vanish in the flat space 
limit, as discussed in Appendix A, as a consequence of stress-energy balance. 

In this generic case eq. (4.19) then implies the following size for the near-source dilaton 
derivative 


lim 


n mV ~^(Ai oc + T) v ~^(T} v ~kV, (5.4) 

>0 L J p v 27 t v Zir v 

with a similar size predicted by (4.26) for the near-brane form of B' (and so also the brane 

defect angle) 

1 - lim (B'W d '\ ~ ~ q ioc) ~ k 2 v 2 . (5.5) 

By contrast, (4.28) predicts the near-brane limit of W’ is additionally suppressed, 

.2 


lim 


¥0 L 


B(W d y — — — (Ai oc > ~ k a v 

j Pv 7 T 


4„,4 


(5.6) 


All of these estimates are confirmed by the numerical results displayed in Fig. 2 and Fig. 6, 
and they are consistent with the weakly gravitating solutions to the Kasner equations (which 
state that 1 — 6 and dw are order z 2 when z<l). 

The d-dimensional curvature in this generic case is not particularly small, since it is of 
order 

r = zi4{y)tot ~ > ( 5 - 7 ) 

corresponding to a vacuum energy of order v 2 . 
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5.2 Scale invariance 

The other extreme is the scale invariant case: (p, q, r) = (—1,1, — 1). As noted previously, in 
this case the quantity y everywhere vanishes, 


y = o, 


(5.8) 


and so the vortex constraints, (4.35) and (4.22), then also imply 


<^loc) v ^ <Aioc) v - 0 . 


(5.9) 


This in turn leads to the a vanishing near-brane limits for both cf>' and W', 


lim 

r v —>0 . 


BW 


d ±t 




and 


lim 

Tv —. 


B{W' 


d\f 


—(A'loc^-O, (5.10) 

- pv 7T 


although the brane defect angle is again given by (5.5) and so is not particularly suppressed. 

This suppression of W' in the near-brane limit resembles the pure Maxwell-Einstein 
(dilaton-free) case considered in [4], for which the radial Einstein constraint also generically 
forces (dfi oc ) v and the near-brane limit of W' to vanish. It is also consistent with the observa¬ 
tion that, in the scale invariant case, there exists the BPS choice of couplings (3 = 1 for which 
X\ oc = Z\ oc = 0 locally and to all orders in r v /£. This is particularly obvious in the numerical 
BPS solution of Fig. 4 which has constant warp factor and dilaton. 

Lastly, the scale-invariant choice ensures the vanishing of the d-dimensional curvature 




(5.11) 


although this is less interesting than it sounds since this is typically achieved by having the 
dilaton zero mode run away, —> 0, since this is the only generic solution to the zero-mode 
equation 0 = {X + T)tot = (A’) to t oc 


Special case: supersymmetric ‘BPS’ branes 

The coefficient of the scale invariant runaway potential vanishes if the flux quantization con¬ 
dition can be satisfied such that Q = Q SUS y, where 

Qsusy := ■ (5.12) 

Kj 

In other words, (A’)tot can be made to vanish in a (^-independent way if this relation holds, 
because it ensures that X B = V B — L A = 0 by having the individual contributions cancel 
against one another locally. The localized contributions to (<T) tot also vanish locally if $ = 1 
and are otherwise suppressed. 

In the scale invariant case, for which <fi' = 0 at the source, we can insert this desired value 
of Q = Qsusy into the the flux quantization condition in (2.45), and derive a local condition 
relating tension and localized flux. For N = +1 it reads 

®a± = — (l-a±) => K 2 T ± = -2g R t ± (5.13) 

9r 
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so long as a+ = a_. Otherwise, the defect angles contribute in a nonlinear way to the flux 
quantization condition via T := ^/a + a_ and no such local condition can be found. Note also 
that the relation (5.12) can only be made to hold if the gauge symmetry has the coupling 
strength of the ii-symmetry, g A = g R . 

Such a relation between the defect angle and flux is also expected from the supersymmetry 
conditions, in order to guarantee the continued cancellation of spin and gauge connections 
within the Killing spinor’s covariant derivative, once localized sources are introduced [5, 17]. 
However, past works did not properly identify the defect angle with the renormalized brane 
tension, as in (5.13). 

5.3 Decoupling case: p = — q = 2r + 1 

The ‘decoupling’ choice p = — q = 2r + 1 is special because it ensures that j3 is ^-independent. 
This suppresses the dilaton-dependence of the brane tension, as in (3.68), which also makes 
it similar in form to the (^-dependence of the brane-localized flux, (. Notice both become 
^-independent in the scale-invariant special case where r = — 1. 

We here ask whether this suppression of brane-dilaton couplings also suppresses the 
brane’s gravitational response, and if so by how much. To decide, recall that for the decoupling 
choice, y takes on the special form 

y = (q — l)(Vh — L z ) + (r + 1 )L blf = (q — l)7fi oc + (r + 1 )L blf , (5-14) 

and the vortex constraint in (4.35) becomes 

« 2 <^ioc) v - k 2 (*Ioc) v - (qX\oc + {r + 1 )L blf ) 2 v ~ (r + 1 ) 2 K 4 (L BiF ) 2 , (5.15) 

which uses that (4.35) requires {X\ oc ) v <C ( x L blf ) v . 

Conveniently, the quantity (L blf } v can be estimated using the on-shell value of the 
4-form field strength (3.44) and the bulk gauge field strength (3.75). This gives 

<W} V = £ Q J d pZpoe^+W ~ -Q e (r+1)0v , (5-16) 

where n v and cj) v are the flux quantum and approximately constant value of the dilaton 
in in the vortex region X v . To get a handle on its size we use the source free estimate 
Q ~ ±g R /n 2 and write the result more transparently in terms of r~ l = e((f>)v = eve~ p ^^ 2 
and = 2 g R (cj))/K = (2 g R /K)e~^/ 2 , to find 

k 2 (L blf ) v ~ -27t n v e e (r+1) ^’ « =F47rn v e = T4vt n v e kv . (5.17) 

For small vortices, r v /£ <C 1, this reveals the decoupling case to lie between the generic and 
scale-invariant cases, with n 2 (y) v suppressed by a single power of the vortex size in KK units, 
as might be expected given that the leading ^-dependence of the point brane action arises 
within the single-derivative localized-flux term of the point-brane action. 
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These estimates lead to the following expectations for near-brane bulk derivatives. As 
always, B' near the sources is dominated by the energy density, with 

„2 


1 — lim 

r v —>-0 


B'W d 


i Pv - 27T^ loc ^’ 


while (by contrast) the near-source derivatives (/)' and W' are KK suppressed 


lim 

r v —>0 - 


BW 


dil 


J Pv 


- - =F( r + !) 2ra v£ ( y ) , 


and 


lim 

r v —>0 . 


R(W d )' 


J Pv 


— (r + l) 2 ( 2n v envY ( 


(5.18) 


(5.19) 


(5.20) 


Because it is suppressed by two powers of f v /£ this last expression is of the same size as terms 
we have neglected, such as second-derivative terms in the brane action, so we should not trust 
its precise numerical prefactor. 

The d-dimensional curvature in this case can be written similarly to give R = 2/c^(y) , 

whose size corresponds to an effective vacuum energy, U e s, of order 


U eS ~ 0>) tot ^ T4tT n v EKV (j) ^ 2 - 


(5.21) 


Strictly speaking, this expression is somewhat self-referential because it is a function of the 
would-be dilaton zero-mode, ip, whose value must be obtained by minimizing a quantity like 
(5.21). Although this generically leads to runaway behaviour in the scale-invariant case (with 
U e ff oc e 2lf implying the minimum occurs for ip —> — oo, the same need not be true when 
r / -1 since the vortex action then breaks scale invariance and so changes the functional 
form of U e ff ((/)). 

What the above expressions leave open is what this precise form is, since this requires 
a more detailed evaluation of the ^-dependence of all quantities that has been done here, 
including all of the ^-dependence implicit within (A) tot , and not just the near-source part 
(A’ioc)v estimated here. Although this takes us beyond the scope of this (already long) paper, 
we do describe such a more detailed calculation in [33], including a description of the 4D 
perspective obtained by integrating out the extra dimensions entirely. 


6. Discussion 

This paper’s aim has been to carefully determine the way in which codimension-two objects 
back-react on their in environment, for a specific UV completion which captures the physics 
of brane-localized flux coupled to the bulk dilaton. To this end, we have determined the way 
in which the microscopic details of the vortex get encoded in IR observables, such as the size 
of the transverse dimensions and the on-brane curvature. 

Quite generally, we find that the breaking of scale invariance in the vortex sector leads 
to modulus stabilization in the IR and — for the particular ‘decoupling’ choice of couplings 
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p = —q = 2r+1 — we find a parametric suppression in the value of the on-brane curvature, by 
a single power of the small ratio r v /i. What remains to be determined is whether reasonable 
choices for vortex-dilaton couplings can stabilize the extra dimensions with a sufficiently large 
hierarchy, i/f v , to profit from this suppression. We explore this in more detail in a companion 
paper [33], where we find such a stabilization to be possible. 

This UV completion verifies earlier claims [6, 8, 10] that moduli can be stabilized in a 
codimension-two version of the Goldberger-Wise mechanism [11], when the dilaton couplings 
of the vortex are chosen to break scale invaria 

This work leaves many open questions. One such asks what the effective description is 
for the dilaton dynamics in the theory below the KK scale within which the extra dimensions 
are integrated out. In particular how does such a theory learn about flux quantization, which 
we’ve seen is central to the dynamics that stabilizes the extra dimensions. We also address 
this question in the companion work [33] . 

Another open direction asks whether vortex configurations can be contrived that break 
supersymmetry in a distributed way (as proposed in [37], for example, with some super- 
symmetry unbroken everywhere locally but with all supersymmetries broken once the entire 
transverse space is taken into account). One might hope to construct a locally half-BPS UV 
vortex — using, eg, a configuration of hyperscalars as in [28] — and embed two (or more) of 
them in the bulk in such a way that leaves supersymmetry completely broken globally. 

We leave these questions for future work. 


NOTE ADDED 

The preprint [38] appeared shortly after our posting of this preprint. As we describe in more 
detail in our reply [39] , we believe its use of 5-function techniques makes it insufficiently precise 
to resolve questions about the size of the vacuum response for 6D systems. It is insufficiently 
precise due to its implicit assumption that the 5-function is independent of the bulk fields; 
an assumption that is most suspicious for the transverse metric, given that the 5-function 
is designed to discriminate points based on their proper distance from the brane sources. 
A proper statement of the bulk field dependence requires a more precise regularization (and 
renormalization) of the brane action, as given for instance in [4, 9] and this paper. The upshot 
of these more precise procedures is that it is stress-energy balance (as expressed through the 
radial Einstein ‘constraint’ equation) that determines any hidden dependence on bulk fields, 
and this is easy to get wrong if one uses a naive regularization for which the regularization 
radius is not set dynamically using a microscopic set of field equations. 
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A. Scaling and the suppression of (X\ oc ) 


We here derive a useful integral identity that is satisfied by the vortex solutions in the limit 
where gravitational back-reaction is neglected so the vortex is in flat space. It is this identity 
that underlies the small size of vortex integrals like (Ti oc ) encountered in the main text. 

The starting point is the observation that the static vortex solution minimizes the energy 
(or negative action) 


I — J V - 5 (^L(p + Lq, + 14 + L a + L z + 


= 2vr 


V^g 


7^2 9 mnd rn(l> d n <f> + V B ((f)) 


+ ^ g mn (d m *d n * + eH 2 Z m Z n ) + j (T 2 - v 2 f 

I I p-<t> A A mn I I pP<t>7 z mn _I_ I PP r <t> 7, 
1 ^ c ' 4 ^ Zjrnn ^ ' 2 


, (A.l) 


and observes that this is stationary with respect to arbitrary variations of the matter fields 
(without also varying the metric), due to their field equations. In particular it is invariant 
under rescalings of the form (j)(y ) -* (f)(y/~s y), d>(y) -> ^(y^ y) and so on. 

Now, suppose the metric g mn also satisfies g m n(y/s y) = sg mn (y), such as is true for the 
case of a locally flat metric, g mn dy m dy n = dp 2 + a 2 p 2 dd 2 under the rescaling p —> y/s p. 
Since I is invariant under arbitrary coordinate transformations, y m —>• £ m (y), in the 2D 
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directions, provided both the matter fields and metric transform, the stationarity of I with 
respect to redefinitions (f>(y) — > (j)(y/s y ) (for all matter fields) is equivalent (for 2D metrics 
with conformal Killing vectors) to stationarity with respect to the rescaling g mn —> sg mn 
without also performing the coordinate rescaling. 

Under the rescaling g mn — » s g mn we have \f--g —» S\[—Q and so (assuming all fields vary 
and point only in the transverse dimensions) 


i — J d yV-ayL'j, + l ^ +14 + l a + l z + L m 

J d 2 y\f~-g L<f> + L\ i> + s (Vb + Vbj + — ^ L A + L z + U m j x ^ 

and so the stationarity condition becomes 

^yV~9 (Vb + ^ L a + L z + L m 




= / d 2 y^ 


U s + V b - (L a + L 


(A.2) 


= / d 2 y^ = 5'*’. 


(A.3) 


The claim is that this equation is an automatic consequence of the matter equations of 
motion, and expresses the balancing of pressures (on average) in the radial directions for a 
stable vortex configuration. 

The same arguments apply equally well for an isolated Q = 0 vortex for which and 
L a are negligible in I, in which case eq. (A.3) reduces to 

0 = J d 2 y^(v b - L z ) = J d 2 y^~g A loc . (A.4) 


In the special BPS case examined in the body this is not only true on average but is also 
locally true, following directly from eq. (3.73). 

Notice that the way it has been derived shows that eq. (A.3) is a statement about the 
vanishing of the extra-dimensional components of the stress energy, as is made more explicit 
in the appendix of [4], It need not hold once the metric back-reaction is turned on, but 
the vanishing of the flat-space result leads to the exact result being is smaller than might 
otherwise have been expected. 
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